i=1
for(i in 1:length(dependentlist)){
dep=dependentlist[i]
print(dependentlist_eng[i])
load(paste("../data/modeldata/RForigmodeldatawithgeoandgeom_",dep,".RData",sep=""))
preds <- evaluateforwardCV_anyerror(mypath=paste("../data/FSCV/RF_geoandgeom/RFwithgeoandgeom_fw_5fold_6p_",dep,"_allpreds",sep=""),kk=1:5,endround = 6,error = "cverror",geheim = "geheimerprederror",yrange=c(0,1))
print(preds)
importance_ranfor_pset(modeldata=origmodeldata,dependent=dep,pset=5,altdata=origmodeldata)
predictors <- c(as.character(preds[1,1]),as.character(preds[2,1]))
predict_ranfor_full(modeldata=origmodeldata, dependent=dep, predictors=predictors,pset=5,kappasum = F,tausum = F,altdata=origmodeldata)
predictors = vector()
for (n in 1:ncol(preds)){
for (n2 in 1:nrow(preds)){
predictors <- c(predictors,as.character(preds[n2,n]))
}
}
print("###################### WITH PREDICTORS from the FW SELECTION ###################")
uniquepredictors <- unique(predictors)
predict_ranfor_full(modeldata=origmodeldata, dependent=dep, predictors=predictors,pset=5,kappasum = F,tausum = F,altdata=origmodeldata)
}
## [1] "1a.2.1 - Potential as a habitat for drought-tolerant species"
## [1] "Prediction error at end is: 0.683982683982684"
## [2] "Prediction error at end is: 0.60995670995671"
## [3] "Prediction error at end is: 0.637662337662338"
## [4] "Prediction error at end is: 0.628571428571429"
## [5] "Prediction error at end is: 0.61948051948052"
## [6] "Prediction error at end is: 0.60995670995671"
## k 1
## 1 geom_10m_fl4_L80
## 2 profc_ws7_hr_hr
## 3 geom_10m_fl8_L8
## 4 geom_hr_L3_fl1_rpliprichness_UE_hr_20cells_hr
## 5 vectorruggedness_hr_ws5
## 6 geom_hr_L50m_fl10_rplippatchdensity_UE_hr_40cells_hr
## k 2
## 1 SGU
## 2 geom_hr_L3_fl1_rpliprichness_UE_hr_20cells_hr
## 3 slope_ws7_10m
## 4 geom_10m_fl3_L22
## 5 geom_10m_fl1_L20
## 6 geom_10m_fl2_L20
## k 3
## 1 geom_10m_fl3_L80
## 2 geom_10m_fl3_L3
## 3 sagaTopographic_Wetness_Index_50m
## 4 dtm_hr_TEXTURE_r10
## 5 crosc_ws3_hr_hr
## 6 geom_hr_L3_fl1_rplipsimpson_UE_hr_20cells_hr
## k 4
## 1 dtm_hr_TEXTURE_r10
## 2 geom_10m_fl3_L50
## 3 geom_hr_L3_fl1_rplippatchdensity_UE_hr_20cells_hr
## 4 geom_10m_fl1_L19
## 5 TRI_hr_ws52
## 6 geom_10m_fl1_L25
## k 5
## 1 geom_10m_fl3_L7
## 2 geom_hr_L50m_fl1_rpliprichness_UE_hr_20cells_hr
## 3 geom_hr_L50m_fl1_rplipshape_UE_hr_40cells_hr
## 4 Profile_Curvature_50m
## 5 geom_hr_L3_fl1_rplipshape_UE_hr_60cells_hr
## 6 maxic_ws11_10m
## [1] "OBB error with all predictors of allpreds is 0.611111111111111"
## MeanDecreaseGini MeanDecreaseAccuracy
## TRI_hr_ws31 0.2823200 0.001872578
## vectorruggedness_hr_ws15 0.3234946 0.001773886
## sagaTopographic_Wetness_Index_50m 0.4709617 0.001348061
## TRI_hr_ws30 0.2515826 0.001338887
## geom_10m_fl4_L9 0.3161175 0.001287507
## Slope_10m 0.1602523 0.001238871
## TRI_hr_ws18 0.2284034 0.001194503
## TRI_hr_ws14 0.2118763 0.001185284
## vectorruggedness_hr_ws59 0.2145666 0.001178964
## TRI_hr_ws35 0.1562579 0.001171159
##
## altpreds 1 2 3 4 5
## 1 4 0 0 0 0
## 2 0 13 0 0 0
## 3 0 0 21 0 0
## 4 0 0 0 38 0
## 5 0 0 0 0 32
## [1] "classification error rate with altdata: 0"
## [1] "OBB error with all predictors of allpreds is 0.583333333333333"
## [1] "OOB-error: 0.527777777777778 for predictors geom_10m_fl4_L80 AND profc_ws7_hr_hr"
##
## preds 1 2 3 4 5
## 1 2 0 0 0 0
## 2 1 9 0 1 2
## 3 0 0 10 0 0
## 4 1 3 11 36 8
## 5 0 1 0 1 22
## [1] "Kappa overall = 0.62114430869723"
## [1] "Tau overall = 0.664351851851852"
## [1] "mean quality = 0.5485649997704"
## [1] "The quality is 0.5485649997704"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 0.6794198163761"
## [1] "classification error rate with altdata: 0.268518518518519"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of allpreds is 0.592592592592593"
## [1] "OOB-error: 0.546296296296296 for predictors geom_10m_fl4_L80 AND profc_ws7_hr_hr AND geom_10m_fl8_L8 AND geom_hr_L3_fl1_rpliprichness_UE_hr_20cells_hr AND vectorruggedness_hr_ws5 AND geom_hr_L50m_fl10_rplippatchdensity_UE_hr_40cells_hr AND SGU AND geom_hr_L3_fl1_rpliprichness_UE_hr_20cells_hr AND slope_ws7_10m AND geom_10m_fl3_L22 AND geom_10m_fl1_L20 AND geom_10m_fl2_L20 AND geom_10m_fl3_L80 AND geom_10m_fl3_L3 AND sagaTopographic_Wetness_Index_50m AND dtm_hr_TEXTURE_r10 AND crosc_ws3_hr_hr AND geom_hr_L3_fl1_rplipsimpson_UE_hr_20cells_hr AND dtm_hr_TEXTURE_r10 AND geom_10m_fl3_L50 AND geom_hr_L3_fl1_rplippatchdensity_UE_hr_20cells_hr AND geom_10m_fl1_L19 AND TRI_hr_ws52 AND geom_10m_fl1_L25 AND geom_10m_fl3_L7 AND geom_hr_L50m_fl1_rpliprichness_UE_hr_20cells_hr AND geom_hr_L50m_fl1_rplipshape_UE_hr_40cells_hr AND Profile_Curvature_50m AND geom_hr_L3_fl1_rplipshape_UE_hr_60cells_hr AND maxic_ws11_10m"
##
## preds 1 2 3 4 5
## 1 4 0 0 0 0
## 2 0 13 0 0 0
## 3 0 0 21 0 0
## 4 0 0 0 38 0
## 5 0 0 0 0 32
## [1] "Kappa overall = 1"
## [1] "Tau overall = 1"
## [1] "mean quality = 1"
## [1] "The quality is 1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 1"
## [1] "classification error rate with altdata: 0"
## [1] "1a.2.2 - Potential as a habitat for moisture-tolerant species"
## [1] "Prediction error at end is: 0.711255411255411"
## [2] "Prediction error at end is: 0.722077922077922"
## [3] "Prediction error at end is: 0.731168831168831"
## [4] "Prediction error at end is: 0.686147186147186"
## [5] "Prediction error at end is: 0.676190476190476"
## [6] "Prediction error at end is: 0.695238095238095"
## k 1
## 1 geom_dtm_10m_hyd_fl5_L6
## 2 vectorruggedness_hr_ws19
## 3 vectorruggedness_hr_ws7
## 4 terraintexture_hr_ws33_tp5
## 5 geom_dtm_10m_hyd_fl5_L37
## 6 geom_10m_fl8_L6
## k 2
## 1 geom_dtm_10m_hyd_fl5_L100
## 2 geom_hr_L3_fl10_rpliprichness_UE_hr_20cells_hr
## 3 SGUcode_vectorruggedness_hr_ws57_TRI_hr_ws31
## 4 minic_ws5_hr_hr
## 5 geom_hr_L50m_fl10_rplipedgedensity_UE_hr_5cells_hr
## 6 geom_10m_fl10_L23
## k 3
## 1 geom_hr_L3_fl1_rplipshannon_UE_hr_20cells_hr
## 2 terraintexture_hr_ws33_tp25
## 3 geom_hr_L3_fl1_rplipsimpson_UE_hr_60cells_hr
## 4 geom_hr_L50m_fl10_rplipshape_UE_hr_5cells_hr
## 5 geom_hr_L3_fl10_rplipshape_UE_hr_5cells_hr
## 6 geom_hr_L50m_fl1_rpliprichness_UE_hr_5cells_hr
## k 4
## 1 geom_10m_fl8_L15
## 2 CrossSectionalCurvature_hr
## 3 Relative_Slope_Position_50m
## 4 geom_10m_fl2_L90
## 5 profc_DTM_50m_avg_ws11_50m
## 6 geom_10m_fl2_L100
## k 5
## 1 dtm_hr_TEXTURE_r10
## 2 terraintexture_hr_ws45_tp25
## 3 geom_hr_L50m_fl1_rplippatchdensity_UE_hr_5cells_hr
## 4 geom_10m_fl4_L14
## 5 geom_hr_L50m_fl1_rplipdominance_UE_hr_10cells_hr
## 6 geom_hr_L50m_fl1_rplipshannon_UE_hr_5cells_hr
## [1] "OBB error with all predictors of allpreds is 0.583333333333333"
## MeanDecreaseGini
## vectorruggedness_hr_ws49 0.1989343
## vectorruggedness_hr_ws27 0.2927074
## Slope_10m 0.1633677
## SGUcode_vectorruggedness_hr_ws57_TRI_hr_ws31.1 0.3949265
## geom_10m_fl4_L9.1 0.2475790
## vectorruggedness_hr_ws25 0.2640938
## TRI_hr_ws21 0.2033363
## terraintexture_hr_ws29_tp25 0.1932838
## maxic_ws7_hr_hr 0.2178423
## vectorruggedness_hr_ws9 0.2809349
## MeanDecreaseAccuracy
## vectorruggedness_hr_ws49 0.001724241
## vectorruggedness_hr_ws27 0.001463012
## Slope_10m 0.001343616
## SGUcode_vectorruggedness_hr_ws57_TRI_hr_ws31.1 0.001274637
## geom_10m_fl4_L9.1 0.001228478
## vectorruggedness_hr_ws25 0.001228461
## TRI_hr_ws21 0.001201919
## terraintexture_hr_ws29_tp25 0.001170891
## maxic_ws7_hr_hr 0.001121009
## vectorruggedness_hr_ws9 0.001067817
##
## altpreds 2 3 4 5
## 2 33 0 0 0
## 3 0 23 0 0
## 4 0 0 35 0
## 5 0 0 0 17
## [1] "classification error rate with altdata: 0"
## [1] "OBB error with all predictors of allpreds is 0.518518518518518"
## [1] "OOB-error: 0.472222222222222 for predictors geom_dtm_10m_hyd_fl5_L6 AND vectorruggedness_hr_ws19"
##
## preds 2 3 4 5
## 2 28 6 9 1
## 3 2 15 1 0
## 4 3 1 25 3
## 5 0 1 0 13
## [1] "Kappa overall = 0.654502369668246"
## [1] "Tau overall = 0.666666666666667"
## [1] "mean quality = 0.616452991452991"
## [1] "The quality is 0.616452991452991"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 0.702116508572032"
## [1] "classification error rate with altdata: 0.25"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of allpreds is 0.537037037037037"
## [1] "OOB-error: 0.462962962962963 for predictors geom_dtm_10m_hyd_fl5_L6 AND vectorruggedness_hr_ws19 AND vectorruggedness_hr_ws7 AND terraintexture_hr_ws33_tp5 AND geom_dtm_10m_hyd_fl5_L37 AND geom_10m_fl8_L6 AND geom_dtm_10m_hyd_fl5_L100 AND geom_hr_L3_fl10_rpliprichness_UE_hr_20cells_hr AND SGUcode_vectorruggedness_hr_ws57_TRI_hr_ws31 AND minic_ws5_hr_hr AND geom_hr_L50m_fl10_rplipedgedensity_UE_hr_5cells_hr AND geom_10m_fl10_L23 AND geom_hr_L3_fl1_rplipshannon_UE_hr_20cells_hr AND terraintexture_hr_ws33_tp25 AND geom_hr_L3_fl1_rplipsimpson_UE_hr_60cells_hr AND geom_hr_L50m_fl10_rplipshape_UE_hr_5cells_hr AND geom_hr_L3_fl10_rplipshape_UE_hr_5cells_hr AND geom_hr_L50m_fl1_rpliprichness_UE_hr_5cells_hr AND geom_10m_fl8_L15 AND CrossSectionalCurvature_hr AND Relative_Slope_Position_50m AND geom_10m_fl2_L90 AND profc_DTM_50m_avg_ws11_50m AND geom_10m_fl2_L100 AND dtm_hr_TEXTURE_r10 AND terraintexture_hr_ws45_tp25 AND geom_hr_L50m_fl1_rplippatchdensity_UE_hr_5cells_hr AND geom_10m_fl4_L14 AND geom_hr_L50m_fl1_rplipdominance_UE_hr_10cells_hr AND geom_hr_L50m_fl1_rplipshannon_UE_hr_5cells_hr"
##
## preds 2 3 4 5
## 2 33 0 0 0
## 3 0 23 0 0
## 4 0 0 35 0
## 5 0 0 0 17
## [1] "Kappa overall = 1"
## [1] "Tau overall = 1"
## [1] "mean quality = 1"
## [1] "The quality is 1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 1"
## [1] "classification error rate with altdata: 0"
## [1] "1a.3 - Habitat for soil organisms"
## [1] "Prediction error at end is: 0.462337662337662"
## [2] "Prediction error at end is: 0.407359307359307"
## [3] "Prediction error at end is: 0.397835497835498"
## [4] "Prediction error at end is: 0.435497835497835"
## [5] "Prediction error at end is: 0.426406926406926"
## [6] "Prediction error at end is: 0.37012987012987"
## k 1
## 1 TRI_hr_ws14
## 2 NormalizedHeight_10m
## 3 slope_ws11_hr_hr
## 4 TRI_hr_ws5
## 5 geom_hr_L50m_fl10_rpliprichness_UE_hr_40cells_hr
## 6 vectorstrength_hr_ws11_hr
## k 2
## 1 SGU
## 2 profc_DTM_50m_avg_ws3_50m
## 3 geom_10m_fl2_L5
## 4 dtm_hr_TEXTURE_r10
## 5 SGUcode_vectorruggedness_hr_ws57_TRI_hr_ws31
## 6 geom_hr_L3_fl1_rplipdominance_UE_hr_40cells_hr
## k 3
## 1 SGU
## 2 geom_10m_fl3_L33
## 3 slope_DTM_50m_avg_ws5_50m
## 4 Channel_Network_Base_Level_hr
## 5 geom_dtm_10m_hyd_fl5_L5
## 6 geom_10m_fl3_L6
## k 4
## 1 SGU
## 2 minic_DTM_50m_avg_ws9_50m
## 3 geom_10m_fl4_L40
## 4 geom_10m_fl3_L33
## 5 vectorruggedness_hr_ws23
## 6 geom_hr_L3_fl1_rplipdominance_UE_hr_40cells_hr
## k 5
## 1 SGU
## 2 elev_ws15_10m
## 3 geom_hr_L50m_fl10_rplippatchdensity_UE_hr_10cells_hr
## 4 geom_hr_L50m_fl10_rplipmps_UE_hr_20cells_hr
## 5 geom_hr_L3_fl1_rpliprichness_UE_hr_5cells_hr
## 6 terraintexture_hr_ws29_tp5
## [1] "OBB error with all predictors of allpreds is 0.37962962962963"
## MeanDecreaseGini MeanDecreaseAccuracy
## elev_ws15_10m 0.8549385 0.006913781
## elev_ws3_10m 0.9015790 0.005765188
## elev_ws5_10m 0.6681540 0.004935603
## Normalized_Height_50m 0.6191737 0.004833930
## vectorruggedness_hr_ws7 0.6394537 0.004193813
## elev_ws11_10m 0.5455804 0.003863798
## ChannelNetworkBaseLevel_10m 0.4663922 0.003795904
## SGU 0.6421372 0.003768861
## vectorruggedness_hr_ws27 0.5064353 0.003511857
## NormalizedHeight_10m 0.5839938 0.003378291
##
## altpreds 1 2 3 4
## 1 24 0 0 0
## 2 0 34 0 0
## 3 0 0 26 0
## 4 0 0 0 24
## [1] "classification error rate with altdata: 0"
## [1] "OBB error with all predictors of allpreds is 0.388888888888889"
## [1] "OOB-error: 0.435185185185185 for predictors TRI_hr_ws14 AND NormalizedHeight_10m"
##
## preds 1 2 3 4
## 1 24 0 0 0
## 2 0 34 0 0
## 3 0 0 26 0
## 4 0 0 0 24
## [1] "Kappa overall = 1"
## [1] "Tau overall = 1"
## [1] "mean quality = 1"
## [1] "The quality is 1"
## [1] "######### Cramer's V = 1"
## [1] "classification error rate with altdata: 0"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of allpreds is 0.416666666666667"
## [1] "OOB-error: 0.259259259259259 for predictors TRI_hr_ws14 AND NormalizedHeight_10m AND slope_ws11_hr_hr AND TRI_hr_ws5 AND geom_hr_L50m_fl10_rpliprichness_UE_hr_40cells_hr AND vectorstrength_hr_ws11_hr AND SGU AND profc_DTM_50m_avg_ws3_50m AND geom_10m_fl2_L5 AND dtm_hr_TEXTURE_r10 AND SGUcode_vectorruggedness_hr_ws57_TRI_hr_ws31 AND geom_hr_L3_fl1_rplipdominance_UE_hr_40cells_hr AND SGU AND geom_10m_fl3_L33 AND slope_DTM_50m_avg_ws5_50m AND Channel_Network_Base_Level_hr AND geom_dtm_10m_hyd_fl5_L5 AND geom_10m_fl3_L6 AND SGU AND minic_DTM_50m_avg_ws9_50m AND geom_10m_fl4_L40 AND geom_10m_fl3_L33 AND vectorruggedness_hr_ws23 AND geom_hr_L3_fl1_rplipdominance_UE_hr_40cells_hr AND SGU AND elev_ws15_10m AND geom_hr_L50m_fl10_rplippatchdensity_UE_hr_10cells_hr AND geom_hr_L50m_fl10_rplipmps_UE_hr_20cells_hr AND geom_hr_L3_fl1_rpliprichness_UE_hr_5cells_hr AND terraintexture_hr_ws29_tp5"
##
## preds 1 2 3 4
## 1 24 0 0 0
## 2 0 34 0 0
## 3 0 0 26 0
## 4 0 0 0 24
## [1] "Kappa overall = 1"
## [1] "Tau overall = 1"
## [1] "mean quality = 1"
## [1] "The quality is 1"
## [1] "######### Cramer's V = 1"
## [1] "classification error rate with altdata: 0"
## [1] "1a.4 - Habitat for crops"
## [1] "Prediction error at end is: 0.177056277056277"
## [2] "Prediction error at end is: 0.195670995670996"
## [3] "Prediction error at end is: 0.222943722943723"
## [4] "Prediction error at end is: 0.240692640692641"
## [5] "Prediction error at end is: 0.232034632034632"
## [6] "Prediction error at end is: 0.232034632034632"
## k 1
## 1 geom_10m_fl10_L5
## 2 slope_ws3_hr_hr
## 3 geom_10m_fl4_L140
## 4 geom_hr_L3_fl10_rplipshape_UE_hr_5cells_hr
## 5 geom_hr_L3_fl10_rpliprichness_UE_hr_60cells_hr
## 6 geom_hr_L3_fl10_rplipedgedensity_UE_hr_5cells_hr
## k 2 k 3
## 1 geom_10m_fl10_L6 geom_10m_fl10_L6
## 2 fischerk_ws27_hr slope_ws5_hr_hr
## 3 fischerk_ws29_hr crosc_DTM_50m_avg_ws5_50m
## 4 planc_DTM_50m_avg_ws3_50m elev_ws11_10m
## 5 PlanCurvature_10m crosc_ws3_hr_hr
## 6 geom_hr_L50m_fl10_rplipmps_UE_hr_5cells_hr CrossSectionalCurvature_50m
## k 4 k 5
## 1 geom_10m_fl10_L6 slope_ws29_hr_hr
## 2 slope_ws3_hr_hr slope_ws5_hr_hr
## 3 geom_10m_fl8_L23 DiurnalAnisotropicHeating_10m
## 4 geom_10m_fl10_L110 Mid_Slope_Positon_50m
## 5 geom_10m_fl10_L140 minic_ws29_hr_hr
## 6 geom_10m_fl3_L32 TRI_hr_ws8
## [1] "OBB error with all predictors of allpreds is 0.157407407407407"
## MeanDecreaseGini MeanDecreaseAccuracy
## slope_ws3_hr_hr 1.2920548 0.011096913
## Slope_hr 1.2351728 0.010711403
## slope_ws5_hr_hr 1.1428811 0.008930896
## TRI_hr_ws1 1.0350235 0.008450871
## MRVBF_hr 0.8705280 0.007082092
## TRI_hr_ws2 0.7726378 0.006599331
## TRI_hr_ws5 0.8846360 0.005547737
## geom_10m_fl10_L6.1 0.6122312 0.004701244
## slope_ws5_10m 0.6775568 0.004596555
## slope_ws19_hr_hr 0.5462017 0.004592360
##
## altpreds 3 4 5
## 3 8 0 0
## 4 0 62 0
## 5 0 0 38
## [1] "classification error rate with altdata: 0"
## [1] "OBB error with all predictors of allpreds is 0.175925925925926"
## [1] "OOB-error: 0.111111111111111 for predictors geom_10m_fl10_L5 AND slope_ws3_hr_hr"
##
## preds 3 4 5
## 3 5 0 0
## 4 2 61 5
## 5 1 1 33
## [1] "Kappa overall = 0.840078973346495"
## [1] "Tau overall = 0.875"
## [1] "mean quality = 0.778019323671498"
## [1] "The quality is 0.778019323671498"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 0.820806690747015"
## [1] "classification error rate with altdata: 0.0833333333333333"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of allpreds is 0.175925925925926"
## [1] "OOB-error: 0.12962962962963 for predictors geom_10m_fl10_L5 AND slope_ws3_hr_hr AND geom_10m_fl4_L140 AND geom_hr_L3_fl10_rplipshape_UE_hr_5cells_hr AND geom_hr_L3_fl10_rpliprichness_UE_hr_60cells_hr AND geom_hr_L3_fl10_rplipedgedensity_UE_hr_5cells_hr AND geom_10m_fl10_L6 AND fischerk_ws27_hr AND fischerk_ws29_hr AND planc_DTM_50m_avg_ws3_50m AND PlanCurvature_10m AND geom_hr_L50m_fl10_rplipmps_UE_hr_5cells_hr AND geom_10m_fl10_L6 AND slope_ws5_hr_hr AND crosc_DTM_50m_avg_ws5_50m AND elev_ws11_10m AND crosc_ws3_hr_hr AND CrossSectionalCurvature_50m AND geom_10m_fl10_L6 AND slope_ws3_hr_hr AND geom_10m_fl8_L23 AND geom_10m_fl10_L110 AND geom_10m_fl10_L140 AND geom_10m_fl3_L32 AND slope_ws29_hr_hr AND slope_ws5_hr_hr AND DiurnalAnisotropicHeating_10m AND Mid_Slope_Positon_50m AND minic_ws29_hr_hr AND TRI_hr_ws8"
##
## preds 3 4 5
## 3 8 0 0
## 4 0 62 0
## 5 0 0 38
## [1] "Kappa overall = 1"
## [1] "Tau overall = 1"
## [1] "mean quality = 1"
## [1] "The quality is 1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 1"
## [1] "classification error rate with altdata: 0"
## [1] "1c.1 - Average precipitation retention capacity"
## [1] "Prediction error at end is: 0.667099567099567"
## [2] "Prediction error at end is: 0.665367965367965"
## [3] "Prediction error at end is: 0.693506493506494"
## [4] "Prediction error at end is: 0.759307359307359"
## [5] "Prediction error at end is: 0.741125541125541"
## [6] "Prediction error at end is: 0.694372294372294"
## k 1
## 1 geom_10m_fl4_L3
## 2 TRI_hr_ws35
## 3 geom_10m_fl10_L8
## 4 geom_10m_fl10_L14
## 5 geom_10m_fl10_L12
## 6 VerticalDistancetoChannelNetwork_10m
## k 2 k 3
## 1 vectorruggedness_hr_ws55 geom_dtm_10m_hyd_fl5_L6
## 2 Vertical_Distance_to_Channel_Network_50m Channel_Network_Base_Level_hr
## 3 vectorruggedness_hr_ws27 dtm_hr_CONVEX_r5
## 4 elev_ws15_10m geom_dtm_10m_hyd_fl5_L3
## 5 geom_10m_fl10_L6 Convexity_10m
## 6 DiurnalAnisotropicHeating_10m crosc_DTM_50m_avg_ws7_50m
## k 4
## 1 crosc_ws23_hr_hr
## 2 geom_10m_fl2_L4
## 3 terraintexture_hr_ws53_tp5
## 4 slope_ws5_hr_hr
## 5 Plan_Curvature_hr
## 6 minic_ws11_10m
## k 5
## 1 terraintexture_hr_ws37_tp5
## 2 terraintexture_hr_ws41_tp5
## 3 DiurnalAnisotropicHeating_10m
## 4 terraintexture_hr_ws29_tp25
## 5 Relative_Slope_Position_hr
## 6 geom_hr_L3_fl10_rpliprichness_UE_hr_60cells_hr
## [1] "OBB error with all predictors of allpreds is 0.62962962962963"
## MeanDecreaseGini
## MaximalCurvature_10m 0.27979971
## vectorruggedness_hr_ws49 0.24700344
## vectorruggedness_hr_ws59 0.33433980
## geom_hr_L3_fl1_rplipdominance_UE_hr_60cells_hr 0.35498139
## geom_10m_fl3_L100 0.08331568
## vectorruggedness_hr_ws37 0.19041487
## geom_10m_fl4_L140.1 0.17079245
## vectorruggedness_hr_ws57 0.18353024
## TRI_hr_ws47 0.24291640
## vectorruggedness_hr_ws35 0.18113058
## MeanDecreaseAccuracy
## MaximalCurvature_10m 0.0011014796
## vectorruggedness_hr_ws49 0.0010956129
## vectorruggedness_hr_ws59 0.0010255557
## geom_hr_L3_fl1_rplipdominance_UE_hr_60cells_hr 0.0009704689
## geom_10m_fl3_L100 0.0009548548
## vectorruggedness_hr_ws37 0.0009332355
## geom_10m_fl4_L140.1 0.0008394995
## vectorruggedness_hr_ws57 0.0008151927
## TRI_hr_ws47 0.0008085507
## vectorruggedness_hr_ws35 0.0007625530
##
## altpreds 1 2 3 4 5
## 1 26 0 0 0 0
## 2 0 45 0 0 0
## 3 0 0 12 0 0
## 4 0 0 0 18 0
## 5 0 0 0 0 7
## [1] "classification error rate with altdata: 0"
## [1] "OBB error with all predictors of allpreds is 0.62962962962963"
## [1] "OOB-error: 0.555555555555556 for predictors geom_10m_fl4_L3 AND TRI_hr_ws35"
##
## preds 1 2 3 4 5
## 1 15 2 3 2 2
## 2 11 43 9 7 3
## 3 0 0 0 0 0
## 4 0 0 0 9 0
## 5 0 0 0 0 2
## [1] "Kappa overall = 0.444253859348199"
## [1] "Tau overall = 0.548611111111111"
## [1] "mean quality = 0.35752380952381"
## [1] "The quality is 0.35752380952381"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = NaN"
## [1] "classification error rate with altdata: 0.361111111111111"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of allpreds is 0.601851851851852"
## [1] "OOB-error: 0.601851851851852 for predictors geom_10m_fl4_L3 AND TRI_hr_ws35 AND geom_10m_fl10_L8 AND geom_10m_fl10_L14 AND geom_10m_fl10_L12 AND VerticalDistancetoChannelNetwork_10m AND vectorruggedness_hr_ws55 AND Vertical_Distance_to_Channel_Network_50m AND vectorruggedness_hr_ws27 AND elev_ws15_10m AND geom_10m_fl10_L6 AND DiurnalAnisotropicHeating_10m AND geom_dtm_10m_hyd_fl5_L6 AND Channel_Network_Base_Level_hr AND dtm_hr_CONVEX_r5 AND geom_dtm_10m_hyd_fl5_L3 AND Convexity_10m AND crosc_DTM_50m_avg_ws7_50m AND crosc_ws23_hr_hr AND geom_10m_fl2_L4 AND terraintexture_hr_ws53_tp5 AND slope_ws5_hr_hr AND Plan_Curvature_hr AND minic_ws11_10m AND terraintexture_hr_ws37_tp5 AND terraintexture_hr_ws41_tp5 AND DiurnalAnisotropicHeating_10m AND terraintexture_hr_ws29_tp25 AND Relative_Slope_Position_hr AND geom_hr_L3_fl10_rpliprichness_UE_hr_60cells_hr"
##
## preds 1 2 3 4 5
## 1 26 0 0 0 0
## 2 0 45 0 0 0
## 3 0 0 12 0 0
## 4 0 0 0 18 0
## 5 0 0 0 0 7
## [1] "Kappa overall = 1"
## [1] "Tau overall = 1"
## [1] "mean quality = 1"
## [1] "The quality is 1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 1"
## [1] "classification error rate with altdata: 0"
## [1] "1c.1 - Minimum precipitation retention capacity"
## [1] "Prediction error at end is: 0.738528138528138"
## [2] "Prediction error at end is: 0.74025974025974"
## [3] "Prediction error at end is: 0.684848484848485"
## [4] "Prediction error at end is: 0.712554112554112"
## [5] "Prediction error at end is: 0.769264069264069"
## [6] "Prediction error at end is: 0.723376623376623"
## k 1
## 1 geom_hr_L3_fl10_rpliprichness_UE_hr_60cells_hr
## 2 maxic_ws13_hr_hr
## 3 geom_hr_L50m_fl1_rplippatchnum_UE_hr_60cells_hr
## 4 geom_10m_fl10_L17
## 5 geom_10m_fl1_L5
## 6 SGUcode_vectorruggedness_hr_ws57_TRI_hr_ws31
## k 2
## 1 terraintexture_hr_ws53_tp25
## 2 SGUcode_vectorruggedness_hr_ws57_TRI_hr_ws31
## 3 Vertical_Distance_to_Channel_Network_50m
## 4 geom_hr_L50m_fl10_rplipdominance_UE_hr_60cells_hr
## 5 geom_hr_L3_fl10_rpliprichness_UE_hr_5cells_hr
## 6 geom_hr_L3_fl1_rplippatchdensity_UE_hr_40cells_hr
## k 3
## 1 terraintexture_hr_ws57_tp25
## 2 terraintexture_hr_ws53_tp25
## 3 geom_dtm_10m_hyd_fl5_L30
## 4 geom_hr_L3_fl10_rplipsimpson_UE_hr_5cells_hr
## 5 geom_hr_L50m_fl1_rplipmps_UE_hr_10cells_hr
## 6 terraintexture_hr_ws17_tp25
## k 4
## 1 profc_ws7_10m
## 2 vectorstrength_hr_ws17_hr
## 3 terraintexture_hr_ws53_tp25
## 4 geom_hr_L3_fl1_rplipdominance_UE_hr_20cells_hr
## 5 minic_ws15_10m
## 6 terraintexture_hr_ws33_tp25
## k 5
## 1 terraintexture_hr_ws57_tp5
## 2 geom_hr_L3_fl10_rplipshannon_UE_hr_5cells_hr
## 3 geom_10m_fl10_L9
## 4 terraintexture_hr_ws41_tp5
## 5 geom_10m_fl4_L17
## 6 geom_hr_L3_fl10_rpliprichness_UE_hr_5cells_hr
## [1] "OBB error with all predictors of allpreds is 0.675925925925926"
## MeanDecreaseGini MeanDecreaseAccuracy
## vectorruggedness_hr_ws11 0.4077027 0.0023059811
## elev_ws5_10m 0.4122884 0.0011638109
## elev_ws11_10m 0.3731010 0.0011110040
## vectorruggedness_hr_ws15 0.3706686 0.0010547322
## vectorruggedness_hr_ws23 0.2011323 0.0010062862
## terraintexture_hr_ws57_tp25 0.3201212 0.0009265658
## TotalCurvature_10m 0.4572167 0.0009247516
## vectorruggedness_hr_ws17 0.2563827 0.0009114961
## Total_Curvature_50m 0.3161605 0.0008903797
## elev_ws3_10m 0.3370111 0.0008665002
##
## altpreds 1 2 3 4 5
## 1 12 0 0 0 0
## 2 0 26 0 0 0
## 3 0 0 22 0 0
## 4 0 0 0 26 0
## 5 0 0 0 0 22
## [1] "classification error rate with altdata: 0"
## [1] "OBB error with all predictors of allpreds is 0.703703703703704"
## [1] "OOB-error: 0.601851851851852 for predictors geom_hr_L3_fl10_rpliprichness_UE_hr_60cells_hr AND maxic_ws13_hr_hr"
##
## preds 1 2 3 4 5
## 1 10 0 0 0 0
## 2 1 25 0 1 1
## 3 0 0 21 1 0
## 4 0 1 0 23 1
## 5 1 0 1 1 20
## [1] "Kappa overall = 0.894071490845684"
## [1] "Tau overall = 0.895833333333333"
## [1] "mean quality = 0.845974869708003"
## [1] "The quality is 0.845974869708003"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 0.898014007665387"
## [1] "classification error rate with altdata: 0.0833333333333333"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of allpreds is 0.685185185185185"
## [1] "OOB-error: 0.537037037037037 for predictors geom_hr_L3_fl10_rpliprichness_UE_hr_60cells_hr AND maxic_ws13_hr_hr AND geom_hr_L50m_fl1_rplippatchnum_UE_hr_60cells_hr AND geom_10m_fl10_L17 AND geom_10m_fl1_L5 AND SGUcode_vectorruggedness_hr_ws57_TRI_hr_ws31 AND terraintexture_hr_ws53_tp25 AND SGUcode_vectorruggedness_hr_ws57_TRI_hr_ws31 AND Vertical_Distance_to_Channel_Network_50m AND geom_hr_L50m_fl10_rplipdominance_UE_hr_60cells_hr AND geom_hr_L3_fl10_rpliprichness_UE_hr_5cells_hr AND geom_hr_L3_fl1_rplippatchdensity_UE_hr_40cells_hr AND terraintexture_hr_ws57_tp25 AND terraintexture_hr_ws53_tp25 AND geom_dtm_10m_hyd_fl5_L30 AND geom_hr_L3_fl10_rplipsimpson_UE_hr_5cells_hr AND geom_hr_L50m_fl1_rplipmps_UE_hr_10cells_hr AND terraintexture_hr_ws17_tp25 AND profc_ws7_10m AND vectorstrength_hr_ws17_hr AND terraintexture_hr_ws53_tp25 AND geom_hr_L3_fl1_rplipdominance_UE_hr_20cells_hr AND minic_ws15_10m AND terraintexture_hr_ws33_tp25 AND terraintexture_hr_ws57_tp5 AND geom_hr_L3_fl10_rplipshannon_UE_hr_5cells_hr AND geom_10m_fl10_L9 AND terraintexture_hr_ws41_tp5 AND geom_10m_fl4_L17 AND geom_hr_L3_fl10_rpliprichness_UE_hr_5cells_hr"
##
## preds 1 2 3 4 5
## 1 12 0 0 0 0
## 2 0 26 0 0 0
## 3 0 0 22 0 0
## 4 0 0 0 26 0
## 5 0 0 0 0 22
## [1] "Kappa overall = 1"
## [1] "Tau overall = 1"
## [1] "mean quality = 1"
## [1] "The quality is 1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 1"
## [1] "classification error rate with altdata: 0"
## [1] "1c.2 - Retention capacity for heavy precipitation events"
## [1] "Prediction error at end is: 0.333333333333333"
## [2] "Prediction error at end is: 0.314718614718615"
## [3] "Prediction error at end is: 0.296536796536797"
## [4] "Prediction error at end is: 0.315584415584416"
## [5] "Prediction error at end is: 0.287445887445887"
## [6] "Prediction error at end is: 0.296969696969697"
## k 1
## 1 geom_10m_fl8_L16
## 2 Modified_Catchment_Area_50m
## 3 Catchment_slope_50m
## 4 geom_hr_L50m_fl10_rpliprichness_UE_hr_60cells_hr
## 5 Catchment_Area2_50m
## 6 Total_Curvature_50m
## k 2
## 1 geom_10m_fl1_L8
## 2 crosc_DTM_50m_avg_ws5_50m
## 3 maxic_ws7_hr_hr
## 4 longc_ws3_hr_hr
## 5 geom_hr_L3_fl1_rplipmps_UE_hr_60cells_hr
## 6 geom_hr_L3_fl1_rplipdominance_UE_hr_20cells_hr
## k 3
## 1 geom_hr_L3_fl10_rplipmps_UE_hr_5cells_hr
## 2 geom_10m_fl2_L8
## 3 crosc_ws11_10m
## 4 geom_hr_L50m_fl10_rplipdominance_UE_hr_5cells_hr
## 5 Catchment_slope_50m
## 6 SlopeHeight_10m
## k 4 k 5
## 1 geom_10m_fl10_L16 geom_10m_fl3_L11
## 2 geom_hr_L3_fl1_rplipshannon_UE_hr_40cells_hr geom_dtm_10m_hyd_fl5_L130
## 3 Vertical_Distance_to_Channel_Network_50m geom_10m_fl3_L120
## 4 fischerk_ws53_hr geom_10m_fl1_L9
## 5 FlowLineCurvature_10m Texture_50m
## 6 vectorruggedness_hr_ws35 geom_dtm_10m_hyd_fl5_L110
## [1] "OBB error with all predictors of allpreds is 0.287037037037037"
## MeanDecreaseGini
## FlowLineCurvature_10m 0.3883731
## vectorruggedness_hr_ws25 0.1386517
## crosc_ws11_10m 0.2424327
## vectorruggedness_hr_ws21 0.2543887
## minic_ws29_hr_hr 0.3022987
## profc_ws15_hr_hr 0.1493986
## maxic_ws15_hr_hr 0.1506604
## vectorruggedness_hr_ws47 0.2067361
## Vertical_Distance_to_Channel_Network_50m 0.2598812
## profc_ws19_hr_hr 0.1860083
## MeanDecreaseAccuracy
## FlowLineCurvature_10m 0.003774068
## vectorruggedness_hr_ws25 0.002113831
## crosc_ws11_10m 0.002012501
## vectorruggedness_hr_ws21 0.001863489
## minic_ws29_hr_hr 0.001776280
## profc_ws15_hr_hr 0.001477575
## maxic_ws15_hr_hr 0.001384619
## vectorruggedness_hr_ws47 0.001379671
## Vertical_Distance_to_Channel_Network_50m 0.001352361
## profc_ws19_hr_hr 0.001246409
##
## altpreds 1 2 3 4 5
## 1 75 0 0 0 0
## 2 0 11 0 0 0
## 3 0 0 6 0 0
## 4 0 0 0 8 0
## 5 0 0 0 0 8
## [1] "classification error rate with altdata: 0"
## [1] "OBB error with all predictors of allpreds is 0.296296296296296"
## [1] "OOB-error: 0.305555555555556 for predictors geom_10m_fl8_L16 AND Modified_Catchment_Area_50m"
##
## preds 1 2 3 4 5
## 1 75 5 6 2 5
## 2 0 6 0 0 0
## 3 0 0 0 0 0
## 4 0 0 0 6 0
## 5 0 0 0 0 3
## [1] "Kappa overall = 0.572841133816744"
## [1] "Tau overall = 0.791666666666667"
## [1] "mean quality = 0.495381231671554"
## [1] "The quality is 0.495381231671554"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = NaN"
## [1] "classification error rate with altdata: 0.166666666666667"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of allpreds is 0.296296296296296"
## [1] "OOB-error: 0.25 for predictors geom_10m_fl8_L16 AND Modified_Catchment_Area_50m AND Catchment_slope_50m AND geom_hr_L50m_fl10_rpliprichness_UE_hr_60cells_hr AND Catchment_Area2_50m AND Total_Curvature_50m AND geom_10m_fl1_L8 AND crosc_DTM_50m_avg_ws5_50m AND maxic_ws7_hr_hr AND longc_ws3_hr_hr AND geom_hr_L3_fl1_rplipmps_UE_hr_60cells_hr AND geom_hr_L3_fl1_rplipdominance_UE_hr_20cells_hr AND geom_hr_L3_fl10_rplipmps_UE_hr_5cells_hr AND geom_10m_fl2_L8 AND crosc_ws11_10m AND geom_hr_L50m_fl10_rplipdominance_UE_hr_5cells_hr AND Catchment_slope_50m AND SlopeHeight_10m AND geom_10m_fl10_L16 AND geom_hr_L3_fl1_rplipshannon_UE_hr_40cells_hr AND Vertical_Distance_to_Channel_Network_50m AND fischerk_ws53_hr AND FlowLineCurvature_10m AND vectorruggedness_hr_ws35 AND geom_10m_fl3_L11 AND geom_dtm_10m_hyd_fl5_L130 AND geom_10m_fl3_L120 AND geom_10m_fl1_L9 AND Texture_50m AND geom_dtm_10m_hyd_fl5_L110"
##
## preds 1 2 3 4 5
## 1 75 0 0 0 0
## 2 0 11 0 0 0
## 3 0 0 6 0 0
## 4 0 0 0 8 0
## 5 0 0 0 0 8
## [1] "Kappa overall = 1"
## [1] "Tau overall = 1"
## [1] "mean quality = 1"
## [1] "The quality is 1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 1"
## [1] "classification error rate with altdata: 0"
## [1] "1c.3 - groundwater reformation rate"
## [1] "Prediction error at end is: 0.704761904761905"
## [2] "Prediction error at end is: 0.731601731601732"
## [3] "Prediction error at end is: 0.665367965367965"
## [4] "Prediction error at end is: 0.654978354978355"
## [5] "Prediction error at end is: 0.598701298701299"
## [6] "Prediction error at end is: 0.637229437229437"
## k 1
## 1 geom_10m_fl8_L5
## 2 geom_hr_L50m_fl1_rplipsimpson_UE_hr_20cells_hr
## 3 geom_hr_L50m_fl1_rpliprichness_UE_hr_10cells_hr
## 4 geom_10m_fl8_L4
## 5 geom_10m_fl8_L3
## 6 geom_hr_L3_fl10_rpliprichness_UE_hr_10cells_hr
## k 2
## 1 SGU
## 2 SGUcode_vectorruggedness_hr_ws57_TRI_hr_ws31
## 3 slope_DTM_50m_avg_ws7_50m
## 4 SGUT_wTGnew
## 5 maxic_DTM_50m_avg_ws5_50m
## 6 geom_10m_fl10_L6
## k 3
## 1 SGU
## 2 geom_10m_fl8_L50
## 3 CrossSectionalCurvature_hr
## 4 Vertical_Distance_to_Channel_Network_50m
## 5 terraintexture_hr_ws29_tp5
## 6 TotalCurvature_10m
## k 4
## 1 terraintexture_hr_ws49_tp25
## 2 TRI_hr_ws17
## 3 geom_hr_L3_fl10_rplipmps_UE_hr_20cells_hr
## 4 fischerk_ws33_hr
## 5 profc_ws7_10m
## 6 terraintexture_hr_ws49_tp5
## k 5
## 1 terraintexture_hr_ws37_tp25
## 2 geom_hr_L50m_fl10_rpliprichness_UE_hr_10cells_hr
## 3 geom_10m_fl8_L7
## 4 vectorstrength_hr_ws13_hr
## 5 geom_10m_fl2_L16
## 6 geom_10m_fl1_L110
## [1] "OBB error with all predictors of allpreds is 0.62962962962963"
## MeanDecreaseGini MeanDecreaseAccuracy
## vectorruggedness_hr_ws59 0.3488135 0.001858079
## vectorruggedness_hr_ws43 0.3607745 0.001646025
## vectorruggedness_hr_ws35 0.3500364 0.001586086
## vectorruggedness_hr_ws55 0.4023994 0.001418956
## vectorruggedness_hr_ws13 0.3501757 0.001317820
## vectorruggedness_hr_ws27 0.3841066 0.001313411
## LSFactor_10m 0.1499655 0.001205072
## vectorruggedness_hr_ws51 0.2968573 0.001154607
## SGU.1 0.4128421 0.001094913
## TRI_hr_ws5 0.2039083 0.001091161
##
## altpreds 1 2 3 4 5
## 1 14 0 0 0 0
## 2 0 39 0 0 0
## 3 0 0 14 0 0
## 4 0 0 0 19 0
## 5 0 0 0 0 22
## [1] "classification error rate with altdata: 0"
## [1] "OBB error with all predictors of allpreds is 0.62037037037037"
## [1] "OOB-error: 0.509259259259259 for predictors geom_10m_fl8_L5 AND geom_hr_L50m_fl1_rplipsimpson_UE_hr_20cells_hr"
##
## preds 1 2 3 4 5
## 1 6 0 3 0 1
## 2 4 34 5 2 5
## 3 0 0 5 0 0
## 4 0 0 0 15 0
## 5 4 5 1 2 16
## [1] "Kappa overall = 0.598279669882599"
## [1] "Tau overall = 0.62962962962963"
## [1] "mean quality = 0.513743985632531"
## [1] "The quality is 0.513743985632531"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 0.640367987990471"
## [1] "classification error rate with altdata: 0.296296296296296"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of allpreds is 0.675925925925926"
## [1] "OOB-error: 0.546296296296296 for predictors geom_10m_fl8_L5 AND geom_hr_L50m_fl1_rplipsimpson_UE_hr_20cells_hr AND geom_hr_L50m_fl1_rpliprichness_UE_hr_10cells_hr AND geom_10m_fl8_L4 AND geom_10m_fl8_L3 AND geom_hr_L3_fl10_rpliprichness_UE_hr_10cells_hr AND SGU AND SGUcode_vectorruggedness_hr_ws57_TRI_hr_ws31 AND slope_DTM_50m_avg_ws7_50m AND SGUT_wTGnew AND maxic_DTM_50m_avg_ws5_50m AND geom_10m_fl10_L6 AND SGU AND geom_10m_fl8_L50 AND CrossSectionalCurvature_hr AND Vertical_Distance_to_Channel_Network_50m AND terraintexture_hr_ws29_tp5 AND TotalCurvature_10m AND terraintexture_hr_ws49_tp25 AND TRI_hr_ws17 AND geom_hr_L3_fl10_rplipmps_UE_hr_20cells_hr AND fischerk_ws33_hr AND profc_ws7_10m AND terraintexture_hr_ws49_tp5 AND terraintexture_hr_ws37_tp25 AND geom_hr_L50m_fl10_rpliprichness_UE_hr_10cells_hr AND geom_10m_fl8_L7 AND vectorstrength_hr_ws13_hr AND geom_10m_fl2_L16 AND geom_10m_fl1_L110"
##
## preds 1 2 3 4 5
## 1 14 0 0 0 0
## 2 0 39 0 0 0
## 3 0 0 14 0 0
## 4 0 0 0 19 0
## 5 0 0 0 0 22
## [1] "Kappa overall = 1"
## [1] "Tau overall = 1"
## [1] "mean quality = 1"
## [1] "The quality is 1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 1"
## [1] "classification error rate with altdata: 0"
## [1] "1c.4 - Potential for providing nutrients for plants"
## [1] "Prediction error at end is: 0.25974025974026"
## [2] "Prediction error at end is: 0.333766233766234"
## [3] "Prediction error at end is: 0.315584415584416"
## [4] "Prediction error at end is: 0.296536796536797"
## [5] "Prediction error at end is: 0.306493506493506"
## [6] "Prediction error at end is: 0.324675324675325"
## k 1
## 1 SGU
## 2 geom_10m_fl10_L3
## 3 vectorstrength_hr_ws17_hr
## 4 geom_10m_fl2_L5
## 5 geom_hr_L50m_fl1_rplipedgedensity_UE_hr_60cells_hr
## 6 VerticalDistancetoChannelNetwork_10m
## k 2
## 1 SGU
## 2 geom_10m_fl8_L11
## 3 geom_hr_L3_fl10_rplippatchdensity_UE_hr_60cells_hr
## 4 geom_hr_L3_fl1_rplipsimpson_UE_hr_5cells_hr
## 5 geom_10m_fl2_L4
## 6 sagaTopographic_Wetness_Index_50m
## k 3
## 1 SGU
## 2 Longitudinal_Curvature_50m
## 3 geom_10m_fl1_L9
## 4 geom_hr_L3_fl1_rpliprichness_UE_hr_10cells_hr
## 5 geom_hr_L3_fl1_rplipmps_UE_hr_10cells_hr
## 6 geom_10m_fl1_L4
## k 4
## 1 SGU
## 2 ValleyDepth_10m
## 3 VerticalDistancetoChannelNetwork_10m
## 4 Texture_10m
## 5 terraintexture_hr_ws29_tp25
## 6 geom_hr_L50m_fl1_rplippatchdensity_UE_hr_10cells_hr
## k 5
## 1 SGU
## 2 geom_10m_fl3_L12
## 3 Vertical_Distance_to_Channel_Network_hr
## 4 Catchment_slope_50m
## 5 geom_10m_fl1_L30
## 6 geom_10m_fl3_L23
## [1] "OBB error with all predictors of allpreds is 0.314814814814815"
## MeanDecreaseGini MeanDecreaseAccuracy
## SGU 0.6680496 0.004606939
## vectorruggedness_hr_ws25 0.6748379 0.004380943
## SGU.1 0.6701422 0.004022976
## vectorruggedness_hr_ws27 0.4492908 0.003844611
## minic_ws7_hr_hr 0.3908256 0.002601657
## TotalCurvature_10m 0.5374482 0.002522411
## vectorruggedness_hr_ws29 0.5675637 0.002382132
## vectorruggedness_hr_ws31 0.3422939 0.002213987
## elev_ws11_10m 0.3290646 0.002190745
## vectorruggedness_hr_ws13 0.2573725 0.002075946
##
## altpreds 1 3 5
## 1 65 0 0
## 3 0 11 0
## 5 0 0 32
## [1] "classification error rate with altdata: 0"
## [1] "OBB error with all predictors of allpreds is 0.305555555555556"
## [1] "OOB-error: 0.305555555555556 for predictors SGU AND geom_10m_fl10_L3"
##
## preds 1 3 5
## 1 60 4 4
## 3 0 1 0
## 5 5 6 28
## [1] "Kappa overall = 0.657142857142857"
## [1] "Tau overall = 0.736111111111111"
## [1] "mean quality = 0.521329896608648"
## [1] "The quality is 0.521329896608648"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 0.572763122715766"
## [1] "classification error rate with altdata: 0.175925925925926"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of allpreds is 0.296296296296296"
## [1] "OOB-error: 0.231481481481481 for predictors SGU AND geom_10m_fl10_L3 AND vectorstrength_hr_ws17_hr AND geom_10m_fl2_L5 AND geom_hr_L50m_fl1_rplipedgedensity_UE_hr_60cells_hr AND VerticalDistancetoChannelNetwork_10m AND SGU AND geom_10m_fl8_L11 AND geom_hr_L3_fl10_rplippatchdensity_UE_hr_60cells_hr AND geom_hr_L3_fl1_rplipsimpson_UE_hr_5cells_hr AND geom_10m_fl2_L4 AND sagaTopographic_Wetness_Index_50m AND SGU AND Longitudinal_Curvature_50m AND geom_10m_fl1_L9 AND geom_hr_L3_fl1_rpliprichness_UE_hr_10cells_hr AND geom_hr_L3_fl1_rplipmps_UE_hr_10cells_hr AND geom_10m_fl1_L4 AND SGU AND ValleyDepth_10m AND VerticalDistancetoChannelNetwork_10m AND Texture_10m AND terraintexture_hr_ws29_tp25 AND geom_hr_L50m_fl1_rplippatchdensity_UE_hr_10cells_hr AND SGU AND geom_10m_fl3_L12 AND Vertical_Distance_to_Channel_Network_hr AND Catchment_slope_50m AND geom_10m_fl1_L30 AND geom_10m_fl3_L23"
##
## preds 1 3 5
## 1 65 0 0
## 3 0 11 0
## 5 0 0 32
## [1] "Kappa overall = 1"
## [1] "Tau overall = 1"
## [1] "mean quality = 1"
## [1] "The quality is 1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 1"
## [1] "classification error rate with altdata: 0"
## [1] "1c.5 - Potential as a CO2 sink"
## [1] "Prediction error at end is: 0.541991341991342"
## [2] "Prediction error at end is: 0.466666666666667"
## [3] "Prediction error at end is: 0.42987012987013"
## [4] "Prediction error at end is: 0.382251082251082"
## [5] "Prediction error at end is: 0.410822510822511"
## [6] "Prediction error at end is: 0.392640692640693"
## k 1
## 1 slope_ws19_hr_hr
## 2 fischerk_ws19_hr
## 3 vectorstrength_hr_ws27_hr
## 4 geom_hr_L3_fl1_rplipsimpson_UE_hr_20cells_hr
## 5 minic_DTM_50m_avg_ws11_50m
## 6 TRI_hr_ws18
## k 2
## 1 slope_ws19_hr_hr
## 2 maxic_DTM_50m_avg_ws5_50m
## 3 vectorruggedness_hr_ws7
## 4 geom_hr_L3_fl1_rplipshape_UE_hr_40cells_hr
## 5 geom_10m_fl3_L110
## 6 geom_hr_L3_fl10_rplipdominance_UE_hr_5cells_hr
## k 3
## 1 CrossSectionalCurvature_hr
## 2 TotalCurvature_10m
## 3 Topographic_Wetness_Index_50m
## 4 Channel_Network_Base_Level_hr
## 5 geom_hr_L50m_fl1_rplipdominance_UE_hr_20cells_hr
## 6 longc_ws5_hr_hr
## k 4 k 5
## 1 crosc_ws5_hr_hr SGU
## 2 slope_ws23_hr_hr vectorruggedness_hr_ws35
## 3 geom_dtm_10m_hyd_fl5_L29 MassBalanceIndex_10m
## 4 geom_hr_L3_fl1_rplipshannon_UE_hr_20cells_hr crosc_ws3_10m
## 5 CrossSectionalCurvature_hr SGUT_wTGnew
## 6 geom_10m_fl2_L140 geom_10m_fl4_L19
## [1] "OBB error with all predictors of allpreds is 0.333333333333333"
## MeanDecreaseGini MeanDecreaseAccuracy
## elev_ws3_10m 0.9418425 0.007665726
## vectorruggedness_hr_ws7 0.7755776 0.006588824
## elev_ws15_10m 0.7687863 0.006467315
## vectorruggedness_hr_ws9 0.3699884 0.005604114
## vectorruggedness_hr_ws21 0.5762638 0.005282938
## vectorruggedness_hr_ws3 0.8258676 0.005266290
## elev_ws11_10m 0.6622378 0.005152287
## elev_ws5_10m 0.6380886 0.005086209
## vectorruggedness_hr_ws5 0.6236385 0.004842494
## vectorruggedness_hr_ws31 0.7160756 0.004086461
##
## altpreds 1 2 3 4 5
## 1 51 0 0 0 0
## 2 0 1 0 0 0
## 3 0 0 18 0 0
## 4 0 0 0 33 0
## 5 0 0 0 0 5
## [1] "classification error rate with altdata: 0"
## [1] "OBB error with all predictors of allpreds is 0.305555555555556"
## [1] "OOB-error: 0.314814814814815 for predictors slope_ws19_hr_hr AND fischerk_ws19_hr"
##
## preds 1 2 3 4 5
## 1 51 0 0 0 0
## 2 0 1 0 0 0
## 3 0 0 18 0 0
## 4 0 0 0 33 0
## 5 0 0 0 0 5
## [1] "Kappa overall = 1"
## [1] "Tau overall = 1"
## [1] "mean quality = 1"
## [1] "The quality is 1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 1"
## [1] "classification error rate with altdata: 0"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of allpreds is 0.333333333333333"
## [1] "OOB-error: 0.240740740740741 for predictors slope_ws19_hr_hr AND fischerk_ws19_hr AND vectorstrength_hr_ws27_hr AND geom_hr_L3_fl1_rplipsimpson_UE_hr_20cells_hr AND minic_DTM_50m_avg_ws11_50m AND TRI_hr_ws18 AND slope_ws19_hr_hr AND maxic_DTM_50m_avg_ws5_50m AND vectorruggedness_hr_ws7 AND geom_hr_L3_fl1_rplipshape_UE_hr_40cells_hr AND geom_10m_fl3_L110 AND geom_hr_L3_fl10_rplipdominance_UE_hr_5cells_hr AND CrossSectionalCurvature_hr AND TotalCurvature_10m AND Topographic_Wetness_Index_50m AND Channel_Network_Base_Level_hr AND geom_hr_L50m_fl1_rplipdominance_UE_hr_20cells_hr AND longc_ws5_hr_hr AND crosc_ws5_hr_hr AND slope_ws23_hr_hr AND geom_dtm_10m_hyd_fl5_L29 AND geom_hr_L3_fl1_rplipshannon_UE_hr_20cells_hr AND CrossSectionalCurvature_hr AND geom_10m_fl2_L140 AND SGU AND vectorruggedness_hr_ws35 AND MassBalanceIndex_10m AND crosc_ws3_10m AND SGUT_wTGnew AND geom_10m_fl4_L19"
##
## preds 1 2 3 4 5
## 1 51 0 0 0 0
## 2 0 1 0 0 0
## 3 0 0 18 0 0
## 4 0 0 0 33 0
## 5 0 0 0 0 5
## [1] "Kappa overall = 1"
## [1] "Tau overall = 1"
## [1] "mean quality = 1"
## [1] "The quality is 1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 1"
## [1] "classification error rate with altdata: 0"
## [1] "1d.1 - Potential for retention of heavy metals"
## [1] "Prediction error at end is: 0.444588744588745"
## [2] "Prediction error at end is: 0.481385281385281"
## [3] "Prediction error at end is: 0.407792207792208"
## [4] "Prediction error at end is: 0.398268398268398"
## [5] "Prediction error at end is: 0.361471861471861"
## [6] "Prediction error at end is: 0.38961038961039"
## k 1
## 1 SGU
## 2 vectorruggedness_hr_ws45
## 3 Channel_Network_Base_Level_hr
## 4 Catchment_slope_50m
## 5 geom_10m_fl4_L47
## 6 fischerk_ws17_hr
## k 2
## 1 geom_10m_fl8_L60
## 2 Convexity_50m
## 3 geom_10m_fl8_L50
## 4 Relative_Slope_Position_50m
## 5 geom_hr_L50m_fl10_rpliprichness_UE_hr_20cells_hr
## 6 geom_10m_fl4_L14
## k 3
## 1 geom_10m_fl8_L8
## 2 geom_hr_L3_fl1_rplipsimpson_UE_hr_40cells_hr
## 3 geom_dtm_10m_hyd_fl5_L60
## 4 geom_10m_fl4_L23
## 5 geom_hr_L50m_fl1_rplipshape_UE_hr_20cells_hr
## 6 Valley_Depth_hr
## k 4
## 1 geom_10m_fl8_L12
## 2 geom_hr_L50m_fl1_rplipshannon_UE_hr_20cells_hr
## 3 geom_10m_fl3_L25
## 4 profc_DTM_50m_avg_ws7_50m
## 5 geom_10m_fl1_L50
## 6 geom_hr_L3_fl10_rplipshannon_UE_hr_40cells_hr
## k 5
## 1 SGU
## 2 geom_10m_fl1_L42
## 3 TangentialCurvature_10m
## 4 ChannelNetworkBaseLevel_10m
## 5 geom_hr_L3_fl1_rplipmps_UE_hr_10cells_hr
## 6 vectorruggedness_hr_ws59
## [1] "OBB error with all predictors of allpreds is 0.342592592592593"
## MeanDecreaseGini
## vectorruggedness_hr_ws7 0.3681345
## elev_ws5_10m 0.4273979
## vectorruggedness_hr_ws23 0.3559804
## vectorruggedness_hr_ws41 0.5187288
## elev_ws3_10m 0.7682106
## vectorruggedness_hr_ws39 0.5674367
## vectorruggedness_hr_ws29 0.5086838
## Vertical_Distance_to_Channel_Network_50m 0.4177763
## vectorruggedness_hr_ws55 0.4823247
## vectorruggedness_hr_ws25 0.5943780
## MeanDecreaseAccuracy
## vectorruggedness_hr_ws7 0.003925552
## elev_ws5_10m 0.003555039
## vectorruggedness_hr_ws23 0.003483173
## vectorruggedness_hr_ws41 0.003437154
## elev_ws3_10m 0.003407096
## vectorruggedness_hr_ws39 0.003403081
## vectorruggedness_hr_ws29 0.003312758
## Vertical_Distance_to_Channel_Network_50m 0.003246063
## vectorruggedness_hr_ws55 0.003198155
## vectorruggedness_hr_ws25 0.003189255
##
## altpreds 1 2 3 4 5
## 1 42 0 0 0 0
## 2 0 7 0 0 0
## 3 0 0 10 0 0
## 4 0 0 0 7 0
## 5 0 0 0 0 42
## [1] "classification error rate with altdata: 0"
## [1] "OBB error with all predictors of allpreds is 0.333333333333333"
## [1] "OOB-error: 0.333333333333333 for predictors SGU AND vectorruggedness_hr_ws45"
##
## preds 1 2 3 4 5
## 1 41 2 3 2 5
## 2 0 5 0 0 0
## 3 0 0 6 0 0
## 4 0 0 0 3 0
## 5 1 0 1 2 37
## [1] "Kappa overall = 0.772631578947368"
## [1] "Tau overall = 0.814814814814815"
## [1] "mean quality = 0.661292845640672"
## [1] "The quality is 0.661292845640672"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 0.768075773614234"
## [1] "classification error rate with altdata: 0.148148148148148"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of allpreds is 0.342592592592593"
## [1] "OOB-error: 0.314814814814815 for predictors SGU AND vectorruggedness_hr_ws45 AND Channel_Network_Base_Level_hr AND Catchment_slope_50m AND geom_10m_fl4_L47 AND fischerk_ws17_hr AND geom_10m_fl8_L60 AND Convexity_50m AND geom_10m_fl8_L50 AND Relative_Slope_Position_50m AND geom_hr_L50m_fl10_rpliprichness_UE_hr_20cells_hr AND geom_10m_fl4_L14 AND geom_10m_fl8_L8 AND geom_hr_L3_fl1_rplipsimpson_UE_hr_40cells_hr AND geom_dtm_10m_hyd_fl5_L60 AND geom_10m_fl4_L23 AND geom_hr_L50m_fl1_rplipshape_UE_hr_20cells_hr AND Valley_Depth_hr AND geom_10m_fl8_L12 AND geom_hr_L50m_fl1_rplipshannon_UE_hr_20cells_hr AND geom_10m_fl3_L25 AND profc_DTM_50m_avg_ws7_50m AND geom_10m_fl1_L50 AND geom_hr_L3_fl10_rplipshannon_UE_hr_40cells_hr AND SGU AND geom_10m_fl1_L42 AND TangentialCurvature_10m AND ChannelNetworkBaseLevel_10m AND geom_hr_L3_fl1_rplipmps_UE_hr_10cells_hr AND vectorruggedness_hr_ws59"
##
## preds 1 2 3 4 5
## 1 42 0 0 0 0
## 2 0 7 0 0 0
## 3 0 0 10 0 0
## 4 0 0 0 7 0
## 5 0 0 0 0 42
## [1] "Kappa overall = 1"
## [1] "Tau overall = 1"
## [1] "mean quality = 1"
## [1] "The quality is 1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 1"
## [1] "classification error rate with altdata: 0"
## [1] "1d.2 - Potential for transforming organic contaminants"
## [1] "Prediction error at end is: 0.6"
## [2] "Prediction error at end is: 0.572294372294372"
## [3] "Prediction error at end is: 0.654978354978355"
## [4] "Prediction error at end is: 0.617748917748918"
## [5] "Prediction error at end is: 0.627705627705628"
## [6] "Prediction error at end is: 0.637662337662338"
## k 1
## 1 geom_10m_fl10_L21
## 2 geom_dtm_10m_hyd_fl5_L110
## 3 geom_hr_L50m_fl1_rplipshannon_UE_hr_40cells_hr
## 4 minic_DTM_50m_avg_ws11_50m
## 5 geom_10m_fl10_L150
## 6 maxic_DTM_50m_avg_ws11_50m
## k 2 k 3
## 1 geom_dtm_10m_hyd_fl5_L8 geom_10m_fl10_L90
## 2 geom_hr_L3_fl1_rplipsimpson_UE_hr_40cells_hr elev_ws3_10m
## 3 geom_hr_L3_fl1_rpliprichness_UE_hr_20cells_hr planc_DTM_50m_avg_ws3_50m
## 4 geom_hr_L3_fl10_rplipdominance_UE_hr_5cells_hr geom_10m_fl10_L120
## 5 geom_hr_L3_fl10_rpliprichness_UE_hr_5cells_hr maxic_ws3_hr_hr
## 6 geom_hr_L3_fl1_rpliprichness_UE_hr_10cells_hr geom_10m_fl8_L120
## k 4
## 1 geom_10m_fl10_L22
## 2 PlanCurvature_10m
## 3 geom_hr_L50m_fl1_rplipsimpson_UE_hr_60cells_hr
## 4 geom_10m_fl3_L6
## 5 geom_10m_fl10_L17
## 6 geom_10m_fl10_L19
## k 5
## 1 geom_10m_fl10_L22
## 2 elev_ws15_10m
## 3 geom_10m_fl8_L9
## 4 Convexity_50m
## 5 geom_hr_L3_fl10_rpliprichness_UE_hr_10cells_hr
## 6 minic_ws3_10m
## [1] "OBB error with all predictors of allpreds is 0.527777777777778"
## MeanDecreaseGini
## elev_ws5_10m 0.5027974
## Vertical_Distance_to_Channel_Network_50m 0.5260216
## vectorruggedness_hr_ws27 0.6309566
## vectorruggedness_hr_ws25 0.4315166
## vectorruggedness_hr_ws33 0.3605873
## vectorruggedness_hr_ws37 0.3469170
## elev_ws15_10m 0.3954492
## vectorruggedness_hr_ws31 0.3766044
## Channel_Network_Base_Level_hr 0.2919898
## elev_ws11_10m 0.4945420
## MeanDecreaseAccuracy
## elev_ws5_10m 0.003366026
## Vertical_Distance_to_Channel_Network_50m 0.003336541
## vectorruggedness_hr_ws27 0.002601471
## vectorruggedness_hr_ws25 0.002287444
## vectorruggedness_hr_ws33 0.001991645
## vectorruggedness_hr_ws37 0.001954972
## elev_ws15_10m 0.001914333
## vectorruggedness_hr_ws31 0.001840628
## Channel_Network_Base_Level_hr 0.001834953
## elev_ws11_10m 0.001675054
##
## altpreds 1 2 3 4 5
## 1 10 0 0 0 0
## 2 0 22 0 0 0
## 3 0 0 28 0 0
## 4 0 0 0 9 0
## 5 0 0 0 0 39
## [1] "classification error rate with altdata: 0"
## [1] "OBB error with all predictors of allpreds is 0.537037037037037"
## [1] "OOB-error: 0.5 for predictors geom_10m_fl10_L21 AND geom_dtm_10m_hyd_fl5_L110"
##
## preds 1 2 3 4 5
## 1 4 1 0 0 1
## 2 4 15 7 2 4
## 3 2 3 16 0 3
## 4 0 0 0 2 0
## 5 0 3 5 5 31
## [1] "Kappa overall = 0.491405698139864"
## [1] "Tau overall = 0.537037037037037"
## [1] "mean quality = 0.396153846153846"
## [1] "The quality is 0.396153846153846"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 0.507139419578754"
## [1] "classification error rate with altdata: 0.37037037037037"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of allpreds is 0.555555555555556"
## [1] "OOB-error: 0.453703703703704 for predictors geom_10m_fl10_L21 AND geom_dtm_10m_hyd_fl5_L110 AND geom_hr_L50m_fl1_rplipshannon_UE_hr_40cells_hr AND minic_DTM_50m_avg_ws11_50m AND geom_10m_fl10_L150 AND maxic_DTM_50m_avg_ws11_50m AND geom_dtm_10m_hyd_fl5_L8 AND geom_hr_L3_fl1_rplipsimpson_UE_hr_40cells_hr AND geom_hr_L3_fl1_rpliprichness_UE_hr_20cells_hr AND geom_hr_L3_fl10_rplipdominance_UE_hr_5cells_hr AND geom_hr_L3_fl10_rpliprichness_UE_hr_5cells_hr AND geom_hr_L3_fl1_rpliprichness_UE_hr_10cells_hr AND geom_10m_fl10_L90 AND elev_ws3_10m AND planc_DTM_50m_avg_ws3_50m AND geom_10m_fl10_L120 AND maxic_ws3_hr_hr AND geom_10m_fl8_L120 AND geom_10m_fl10_L22 AND PlanCurvature_10m AND geom_hr_L50m_fl1_rplipsimpson_UE_hr_60cells_hr AND geom_10m_fl3_L6 AND geom_10m_fl10_L17 AND geom_10m_fl10_L19 AND geom_10m_fl10_L22 AND elev_ws15_10m AND geom_10m_fl8_L9 AND Convexity_50m AND geom_hr_L3_fl10_rpliprichness_UE_hr_10cells_hr AND minic_ws3_10m"
##
## preds 1 2 3 4 5
## 1 10 0 0 0 0
## 2 0 22 0 0 0
## 3 0 0 28 0 0
## 4 0 0 0 9 0
## 5 0 0 0 0 39
## [1] "Kappa overall = 1"
## [1] "Tau overall = 1"
## [1] "mean quality = 1"
## [1] "The quality is 1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 1"
## [1] "classification error rate with altdata: 0"
## [1] "1d.3 - Potential as filter and buffer for organic contaminants"
## [1] "Prediction error at end is: 0.0826839826839827"
## [2] "Prediction error at end is: 0.0463203463203463"
## [3] "Prediction error at end is: 0.0463203463203463"
## [4] "Prediction error at end is: 0.0372294372294372"
## [5] "Prediction error at end is: 0.0372294372294372"
## [6] "Prediction error at end is: 0.0372294372294372"
## k 1
## 1 Slope_Height_50m
## 2 geom_10m_fl2_L15
## 3 geom_10m_fl2_L16
## 4 geom_10m_fl2_L17
## 5 geom_10m_fl2_L18
## 6 geom_hr_L3_fl10_rplipedgedensity_UE_hr_40cells_hr
## k 2
## 1 geom_hr_L3_fl10_rplipdominance_UE_hr_60cells_hr
## 2 planc_DTM_50m_avg_ws5_50m
## 3 geom_10m_fl2_L19
## 4 geom_10m_fl2_L22
## 5 geom_10m_fl4_L3
## 6 geom_10m_fl8_L17
## k 3
## 1 planc_ws19_hr_hr
## 2 MinimalCurvature_10m
## 3 crosc_ws11_hr_hr
## 4 geom_hr_L50m_fl1_rpliprichness_UE_hr_20cells_hr
## 5 geom_hr_L50m_fl1_rpliprichness_UE_hr_60cells_hr
## 6 fischerk_ws31_hr
## k 4
## 1 maxic_ws7_hr_hr
## 2 geom_hr_L50m_fl10_rpliprichness_UE_hr_10cells_hr
## 3 geom_hr_L50m_fl10_rpliprichness_UE_hr_60cells_hr
## 4 geom_10m_fl10_L3
## 5 minic_ws15_hr_hr
## 6 geom_10m_fl10_L11
## k 5
## 1 geom_hr_L3_fl10_rplipdominance_UE_hr_60cells_hr
## 2 geom_10m_fl2_L19
## 3 geom_10m_fl4_L8
## 4 geom_10m_fl1_L60
## 5 geom_10m_fl2_L60
## 6 minic_ws7_10m
## [1] "OBB error with all predictors of allpreds is 0.037037037037037"
## MeanDecreaseGini
## geom_hr_L50m_fl1_rplipsimpson_UE_hr_20cells_hr 0.046223453
## geom_10m_fl2_L17.1 0.082940315
## geom_hr_L50m_fl1_rplipshape_UE_hr_20cells_hr 0.029509434
## PlanCurvature_10m 0.054469869
## Mid_Slope_Positon_50m 0.085348233
## geom_10m_fl3_L130.1 0.035999554
## geom_hr_L50m_fl1_rplipmps_UE_hr_40cells_hr 0.071236465
## geom_10m_fl8_L40 0.003010885
## dtm_hr_TEXTURE_r20 0.007466083
## geom_hr_L3_fl1_rplippatchdensity_UE_hr_10cells_hr 0.006166667
## MeanDecreaseAccuracy
## geom_hr_L50m_fl1_rplipsimpson_UE_hr_20cells_hr 0.0004327331
## geom_10m_fl2_L17.1 0.0004074340
## geom_hr_L50m_fl1_rplipshape_UE_hr_20cells_hr 0.0003541687
## PlanCurvature_10m 0.0003217833
## Mid_Slope_Positon_50m 0.0003035496
## geom_10m_fl3_L130.1 0.0002844969
## geom_hr_L50m_fl1_rplipmps_UE_hr_40cells_hr 0.0002628852
## geom_10m_fl8_L40 0.0002554463
## dtm_hr_TEXTURE_r20 0.0002481203
## geom_hr_L3_fl1_rplippatchdensity_UE_hr_10cells_hr 0.0002417582
##
## altpreds 4 5
## 4 4 0
## 5 0 104
## [1] "classification error rate with altdata: 0"
## [1] "OBB error with all predictors of allpreds is 0.037037037037037"
## [1] "OOB-error: 0.0185185185185185 for predictors Slope_Height_50m AND geom_10m_fl2_L15"
##
## preds 4 5
## 4 2 0
## 5 2 104
## [1] "Kappa overall = 0.658227848101266"
## [1] "Tau overall = 0.962962962962963"
## [1] "mean quality = 0.740566037735849"
## [1] "The quality is 0.740566037735849"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 0.518568491248852"
## [1] "classification error rate with altdata: 0.0185185185185185"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of allpreds is 0.037037037037037"
## [1] "OOB-error: 0.0277777777777778 for predictors Slope_Height_50m AND geom_10m_fl2_L15 AND geom_10m_fl2_L16 AND geom_10m_fl2_L17 AND geom_10m_fl2_L18 AND geom_hr_L3_fl10_rplipedgedensity_UE_hr_40cells_hr AND geom_hr_L3_fl10_rplipdominance_UE_hr_60cells_hr AND planc_DTM_50m_avg_ws5_50m AND geom_10m_fl2_L19 AND geom_10m_fl2_L22 AND geom_10m_fl4_L3 AND geom_10m_fl8_L17 AND planc_ws19_hr_hr AND MinimalCurvature_10m AND crosc_ws11_hr_hr AND geom_hr_L50m_fl1_rpliprichness_UE_hr_20cells_hr AND geom_hr_L50m_fl1_rpliprichness_UE_hr_60cells_hr AND fischerk_ws31_hr AND maxic_ws7_hr_hr AND geom_hr_L50m_fl10_rpliprichness_UE_hr_10cells_hr AND geom_hr_L50m_fl10_rpliprichness_UE_hr_60cells_hr AND geom_10m_fl10_L3 AND minic_ws15_hr_hr AND geom_10m_fl10_L11 AND geom_hr_L3_fl10_rplipdominance_UE_hr_60cells_hr AND geom_10m_fl2_L19 AND geom_10m_fl4_L8 AND geom_10m_fl1_L60 AND geom_10m_fl2_L60 AND minic_ws7_10m"
##
## preds 4 5
## 4 4 0
## 5 0 104
## [1] "Kappa overall = 1"
## [1] "Tau overall = 1"
## [1] "mean quality = 1"
## [1] "The quality is 1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 0.870192307692308"
## [1] "classification error rate with altdata: 0"
## [1] "1d.4 - Potential for retention of water-soluble contaminants"
## [1] "Prediction error at end is: 0.583549783549784"
## [2] "Prediction error at end is: 0.582683982683983"
## [3] "Prediction error at end is: 0.601298701298701"
## [4] "Prediction error at end is: 0.583116883116883"
## [5] "Prediction error at end is: 0.582683982683983"
## [6] "Prediction error at end is: 0.564069264069264"
## k 1
## 1 geom_10m_fl8_L5
## 2 Relative_Slope_Position_50m
## 3 minic_DTM_50m_avg_ws5_50m
## 4 geom_10m_fl2_L15
## 5 geom_10m_fl1_L7
## 6 Plan_Curvature_hr
## k 2
## 1 geom_10m_fl3_L120
## 2 geom_hr_L50m_fl1_rplipshannon_UE_hr_20cells_hr
## 3 geom_hr_L50m_fl10_rplipshannon_UE_hr_5cells_hr
## 4 geom_10m_fl3_L130
## 5 geom_hr_L50m_fl10_rplipmps_UE_hr_5cells_hr
## 6 geom_hr_L50m_fl1_rplipedgedensity_UE_hr_5cells_hr
## k 3 k 4
## 1 geom_10m_fl8_L4 geom_10m_fl8_L110
## 2 Relative_Slope_Position_50m longc_DTM_50m_avg_ws11_50m
## 3 TRI_hr_ws53 vectorruggedness_hr_ws35
## 4 maxic_ws5_hr_hr slope_ws15_10m
## 5 geom_10m_fl8_L5 SGUcode_vectorruggedness_hr_ws57_TRI_hr_ws31
## 6 NormalizedHeight_10m profc_ws5_hr_hr
## k 5
## 1 geom_10m_fl10_L14
## 2 geom_10m_fl10_L10
## 3 geom_hr_L3_fl1_rpliprichness_UE_hr_40cells_hr
## 4 geom_hr_L50m_fl10_rplipdominance_UE_hr_5cells_hr
## 5 geom_dtm_10m_hyd_fl5_L32
## 6 geom_hr_L3_fl1_rplippatchdensity_UE_hr_20cells_hr
## [1] "OBB error with all predictors of allpreds is 0.5"
## MeanDecreaseGini MeanDecreaseAccuracy
## sagaTopographic_Wetness_Index_50m 0.4100743 0.002946161
## vectorruggedness_hr_ws35 0.4577494 0.002292830
## vectorruggedness_hr_ws41 0.4793662 0.002267833
## SGU.1 0.2898393 0.002259740
## TRI_hr_ws28 0.2288264 0.001755498
## elev_ws11_10m 0.4273840 0.001746943
## vectorruggedness_hr_ws25 0.3132226 0.001589001
## vectorruggedness_hr_ws23 0.3064092 0.001528978
## vectorruggedness_hr_ws27 0.3221944 0.001405396
## TWI_10m 0.1942418 0.001397361
##
## altpreds 1 2 3 4 5
## 1 47 0 0 0 0
## 2 0 14 0 0 0
## 3 0 0 17 0 0
## 4 0 0 0 22 0
## 5 0 0 0 0 8
## [1] "classification error rate with altdata: 0"
## [1] "OBB error with all predictors of allpreds is 0.481481481481481"
## [1] "OOB-error: 0.425925925925926 for predictors geom_10m_fl8_L5 AND Relative_Slope_Position_50m"
##
## preds 1 2 3 4 5
## 1 46 4 4 4 0
## 2 0 7 1 1 0
## 3 0 0 7 0 0
## 4 1 3 5 17 3
## 5 0 0 0 0 5
## [1] "Kappa overall = 0.649656893325016"
## [1] "Tau overall = 0.699074074074074"
## [1] "mean quality = 0.550785144566301"
## [1] "The quality is 0.550785144566301"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 0.675892376003383"
## [1] "classification error rate with altdata: 0.240740740740741"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of allpreds is 0.518518518518518"
## [1] "OOB-error: 0.5 for predictors geom_10m_fl8_L5 AND Relative_Slope_Position_50m AND minic_DTM_50m_avg_ws5_50m AND geom_10m_fl2_L15 AND geom_10m_fl1_L7 AND Plan_Curvature_hr AND geom_10m_fl3_L120 AND geom_hr_L50m_fl1_rplipshannon_UE_hr_20cells_hr AND geom_hr_L50m_fl10_rplipshannon_UE_hr_5cells_hr AND geom_10m_fl3_L130 AND geom_hr_L50m_fl10_rplipmps_UE_hr_5cells_hr AND geom_hr_L50m_fl1_rplipedgedensity_UE_hr_5cells_hr AND geom_10m_fl8_L4 AND Relative_Slope_Position_50m AND TRI_hr_ws53 AND maxic_ws5_hr_hr AND geom_10m_fl8_L5 AND NormalizedHeight_10m AND geom_10m_fl8_L110 AND longc_DTM_50m_avg_ws11_50m AND vectorruggedness_hr_ws35 AND slope_ws15_10m AND SGUcode_vectorruggedness_hr_ws57_TRI_hr_ws31 AND profc_ws5_hr_hr AND geom_10m_fl10_L14 AND geom_10m_fl10_L10 AND geom_hr_L3_fl1_rpliprichness_UE_hr_40cells_hr AND geom_hr_L50m_fl10_rplipdominance_UE_hr_5cells_hr AND geom_dtm_10m_hyd_fl5_L32 AND geom_hr_L3_fl1_rplippatchdensity_UE_hr_20cells_hr"
##
## preds 1 2 3 4 5
## 1 47 0 0 0 0
## 2 0 14 0 0 0
## 3 0 0 17 0 0
## 4 0 0 0 22 0
## 5 0 0 0 0 8
## [1] "Kappa overall = 1"
## [1] "Tau overall = 1"
## [1] "mean quality = 1"
## [1] "The quality is 1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 1"
## [1] "classification error rate with altdata: 0"
## [1] "1d.5 - Potential as buffer for acidic contaminants"
## [1] "Prediction error at end is: 0.676190476190476"
## [2] "Prediction error at end is: 0.722077922077922"
## [3] "Prediction error at end is: 0.712987012987013"
## [4] "Prediction error at end is: 0.703030303030303"
## [5] "Prediction error at end is: 0.667099567099567"
## [6] "Prediction error at end is: 0.603030303030303"
## k 1 k 2
## 1 geom_dtm_10m_hyd_fl5_L110 geom_dtm_10m_hyd_fl5_L140
## 2 Catchment_area_50m planc_DTM_50m_avg_ws9_50m
## 3 geom_10m_fl4_L40 SGU
## 4 planc_ws5_hr_hr crosc_ws5_hr_hr
## 5 geom_dtm_10m_hyd_fl5_L42 geom_10m_fl10_L3
## 6 terraintexture_hr_ws53_tp25 geom_dtm_10m_hyd_fl5_L4
## k 3
## 1 geom_hr_L3_fl10_rplipshannon_UE_hr_10cells_hr
## 2 geom_dtm_10m_hyd_fl5_L110
## 3 maxic_ws3_hr_hr
## 4 terraintexture_hr_ws41_tp5
## 5 planc_ws29_hr_hr
## 6 LS_Factor_50m
## k 4 k 5
## 1 geom_10m_fl1_L12 geom_10m_fl1_L8
## 2 geom_hr_L3_fl10_rplippatchdensity_UE_hr_20cells_hr geom_10m_fl4_L49
## 3 StandardizedHeight_10m elev_ws5_10m
## 4 geom_dtm_10m_hyd_fl5_L5 Catchment_slope_50m
## 5 geom_10m_fl1_L4 slope_ws29_hr_hr
## 6 geom_hr_L50m_fl10_rpliprichness_UE_hr_5cells_hr FlowLineCurvature_10m
## [1] "OBB error with all predictors of allpreds is 0.518518518518518"
## MeanDecreaseGini
## Normalized_Height_50m 0.8388551
## Modified_Catchment_Area_50m 0.6094262
## Standardized_Height_50m 0.6575744
## Vertical_Distance_to_Channel_Network_50m 0.3790405
## Catchment_Area2_50m 0.3992647
## SGUT_wTGnew 0.4280335
## SGUcode_vectorruggedness_hr_ws57_TRI_hr_ws31.1 0.2957347
## vectorstrength_hr_ws33_hr 0.1977812
## Topographic_Wetness_Index_50m 0.2297360
## SGUT_wTGnew.1 0.4821263
## MeanDecreaseAccuracy
## Normalized_Height_50m 0.005078172
## Modified_Catchment_Area_50m 0.002624196
## Standardized_Height_50m 0.002144330
## Vertical_Distance_to_Channel_Network_50m 0.001967351
## Catchment_Area2_50m 0.001936848
## SGUT_wTGnew 0.001698795
## SGUcode_vectorruggedness_hr_ws57_TRI_hr_ws31.1 0.001624830
## vectorstrength_hr_ws33_hr 0.001556860
## Topographic_Wetness_Index_50m 0.001488524
## SGUT_wTGnew.1 0.001394048
##
## altpreds 1 2 3 4 5
## 1 39 0 0 0 0
## 2 0 24 0 0 0
## 3 0 0 20 0 0
## 4 0 0 0 16 0
## 5 0 0 0 0 9
## [1] "classification error rate with altdata: 0"
## [1] "OBB error with all predictors of allpreds is 0.546296296296296"
## [1] "OOB-error: 0.537037037037037 for predictors geom_dtm_10m_hyd_fl5_L110 AND Catchment_area_50m"
##
## preds 1 2 3 4 5
## 1 37 6 4 0 1
## 2 0 15 1 2 0
## 3 2 3 14 1 3
## 4 0 0 1 13 1
## 5 0 0 0 0 4
## [1] "Kappa overall = 0.686920222634508"
## [1] "Tau overall = 0.710648148148148"
## [1] "mean quality = 0.588996168582375"
## [1] "The quality is 0.588996168582375"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 0.702038641647599"
## [1] "classification error rate with altdata: 0.231481481481481"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of allpreds is 0.583333333333333"
## [1] "OOB-error: 0.416666666666667 for predictors geom_dtm_10m_hyd_fl5_L110 AND Catchment_area_50m AND geom_10m_fl4_L40 AND planc_ws5_hr_hr AND geom_dtm_10m_hyd_fl5_L42 AND terraintexture_hr_ws53_tp25 AND geom_dtm_10m_hyd_fl5_L140 AND planc_DTM_50m_avg_ws9_50m AND SGU AND crosc_ws5_hr_hr AND geom_10m_fl10_L3 AND geom_dtm_10m_hyd_fl5_L4 AND geom_hr_L3_fl10_rplipshannon_UE_hr_10cells_hr AND geom_dtm_10m_hyd_fl5_L110 AND maxic_ws3_hr_hr AND terraintexture_hr_ws41_tp5 AND planc_ws29_hr_hr AND LS_Factor_50m AND geom_10m_fl1_L12 AND geom_hr_L3_fl10_rplippatchdensity_UE_hr_20cells_hr AND StandardizedHeight_10m AND geom_dtm_10m_hyd_fl5_L5 AND geom_10m_fl1_L4 AND geom_hr_L50m_fl10_rpliprichness_UE_hr_5cells_hr AND geom_10m_fl1_L8 AND geom_10m_fl4_L49 AND elev_ws5_10m AND Catchment_slope_50m AND slope_ws29_hr_hr AND FlowLineCurvature_10m"
##
## preds 1 2 3 4 5
## 1 39 0 0 0 0
## 2 0 24 0 0 0
## 3 0 0 20 0 0
## 4 0 0 0 16 0
## 5 0 0 0 0 9
## [1] "Kappa overall = 1"
## [1] "Tau overall = 1"
## [1] "mean quality = 1"
## [1] "The quality is 1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 1"
## [1] "classification error rate with altdata: 0"
i=1
for(i in 1:length(dependentlist)){
dep=dependentlist[i]
print(dependentlist_eng[i])
load(paste("../data/modeldata/RForigmodeldatawithgeoandgeom_",dep,".RData",sep=""))
preds <- evaluateforwardCV_anyerror(mypath=paste("../data/FSCV/RF_localterrain/RFwithgeoandgeom_fw_5fold_6p_",dep,"_localterrain",sep=""),kk=1:5,endround = 6,error = "cverror",geheim = "geheimerprederror",yrange=c(0,1))
print(preds)
importance_ranfor_pset(modeldata=origmodeldata,dependent=dep,pset=1,altdata=origmodeldata,withalt=F)
predictors <- c(as.character(preds[1,1]),as.character(preds[2,1]))
predict_ranfor_full(modeldata=origmodeldata, dependent=dep, predictors=predictors,pset=1,kappasum = F,tausum = F,altdata=origmodeldata)
predictors = vector()
for (n in 1:ncol(preds)){
for (n2 in 1:nrow(preds)){
predictors <- c(predictors,as.character(preds[n2,n]))
}
}
print("###################### WITH PREDICTORS from the FW SELECTION ###################")
uniquepredictors <- unique(predictors)
predict_ranfor_full(modeldata=origmodeldata, dependent=dep, predictors=predictors,pset=1,kappasum = F,tausum = F,altdata=origmodeldata)
}
## [1] "1a.2.1 - Potential as a habitat for drought-tolerant species"
## [1] "Prediction error at end is: 0.703463203463203"
## [2] "Prediction error at end is: 0.731168831168831"
## [3] "Prediction error at end is: 0.704329004329004"
## [4] "Prediction error at end is: 0.686580086580087"
## [5] "Prediction error at end is: 0.667532467532468"
## [6] "Prediction error at end is: 0.685714285714286"
## k 1 k 2
## 1 Slope_50m minic_ws23_hr_hr
## 2 CrossSectionalCurvature_hr profc_ws13_hr_hr
## 3 longc_DTM_50m_avg_ws3_50m slope_DTM_50m_avg_ws11_50m
## 4 MaximalCurvature_10m Convexity_50m
## 5 longc_DTM_50m_avg_ws5_50m minic_ws19_hr_hr
## 6 profc_DTM_50m_avg_ws3_50m profc_ws11_hr_hr
## k 3 k 4
## 1 maxic_ws9_hr_hr crosc_ws3_10m
## 2 crosc_ws7_10m Total_Curvature_50m
## 3 slope_DTM_50m_avg_ws11_50m Plan_Curvature_hr
## 4 planc_ws5_10m planc_DTM_50m_avg_ws3_50m
## 5 maxic_ws5_hr_hr planc_ws3_hr_hr
## 6 minic_ws11_hr_hr Slope_50m
## k 5
## 1 crosc_ws19_hr_hr
## 2 Plan_Curvature_hr
## 3 longc_ws5_hr_hr
## 4 crosc_ws23_hr_hr
## 5 maxic_ws3_hr_hr
## 6 Maximal_Curvature_hr
## [1] "OBB error with all predictors of localterrain is 0.601851851851852"
## MeanDecreaseGini MeanDecreaseAccuracy
## slope_ws3_10m 1.0265466 0.003632591
## slope_ws5_10m 0.9122791 0.003233536
## slope_ws13_hr_hr 0.8457388 0.003166309
## Slope_50m 0.9580178 0.002916099
## slope_ws15_10m 0.7347631 0.002781884
## slope_DTM_50m_avg_ws9_50m 0.8000160 0.002553896
## slope_ws11_hr_hr 0.8156472 0.002542512
## crosc_ws19_hr_hr 0.5710928 0.002400407
## Slope_10m 0.8598335 0.002237079
## maxic_ws11_hr_hr 0.5930079 0.002120002
## [1] "OBB error with all predictors of localterrain is 0.62037037037037"
## [1] "OOB-error: 0.564814814814815 for predictors Slope_50m AND CrossSectionalCurvature_hr"
##
## preds 1 2 3 4 5
## 1 4 0 0 0 0
## 2 0 13 0 0 0
## 3 0 0 21 0 0
## 4 0 0 0 38 0
## 5 0 0 0 0 32
## [1] "Kappa overall = 1"
## [1] "Tau overall = 1"
## [1] "mean quality = 1"
## [1] "The quality is 1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 1"
## [1] "classification error rate with altdata: 0"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of localterrain is 0.611111111111111"
## [1] "OOB-error: 0.574074074074074 for predictors Slope_50m AND CrossSectionalCurvature_hr AND longc_DTM_50m_avg_ws3_50m AND MaximalCurvature_10m AND longc_DTM_50m_avg_ws5_50m AND profc_DTM_50m_avg_ws3_50m AND minic_ws23_hr_hr AND profc_ws13_hr_hr AND slope_DTM_50m_avg_ws11_50m AND Convexity_50m AND minic_ws19_hr_hr AND profc_ws11_hr_hr AND maxic_ws9_hr_hr AND crosc_ws7_10m AND slope_DTM_50m_avg_ws11_50m AND planc_ws5_10m AND maxic_ws5_hr_hr AND minic_ws11_hr_hr AND crosc_ws3_10m AND Total_Curvature_50m AND Plan_Curvature_hr AND planc_DTM_50m_avg_ws3_50m AND planc_ws3_hr_hr AND Slope_50m AND crosc_ws19_hr_hr AND Plan_Curvature_hr AND longc_ws5_hr_hr AND crosc_ws23_hr_hr AND maxic_ws3_hr_hr AND Maximal_Curvature_hr"
##
## preds 1 2 3 4 5
## 1 4 0 0 0 0
## 2 0 13 0 0 0
## 3 0 0 21 0 0
## 4 0 0 0 38 0
## 5 0 0 0 0 32
## [1] "Kappa overall = 1"
## [1] "Tau overall = 1"
## [1] "mean quality = 1"
## [1] "The quality is 1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 1"
## [1] "classification error rate with altdata: 0"
## [1] "1a.2.2 - Potential as a habitat for moisture-tolerant species"
## [1] "Prediction error at end is: 0.769264069264069"
## [2] "Prediction error at end is: 0.546320346320346"
## [3] "Prediction error at end is: 0.564502164502165"
## [4] "Prediction error at end is: 0.584415584415584"
## [5] "Prediction error at end is: 0.584415584415584"
## [6] "Prediction error at end is: 0.583982683982684"
## k 1 k 2
## 1 Slope_50m Plan_Curvature_hr
## 2 GeneralCurvature_10m crosc_ws13_hr_hr
## 3 Flow_Line_Curvature_hr maxic_ws13_hr_hr
## 4 LongitudinalCurvature_10m slope_DTM_50m_avg_ws3_50m
## 5 Minimal_Curvature_hr slope_ws29_hr_hr
## 6 longc_ws5_10m Profile_Curvature_50m
## k 3 k 4
## 1 Total_Curvature_50m maxic_ws9_hr_hr
## 2 Convexity_50m slope_DTM_50m_avg_ws5_50m
## 3 maxic_ws11_10m maxic_ws5_hr_hr
## 4 slope_DTM_50m_avg_ws3_50m maxic_ws19_hr_hr
## 5 profc_ws15_10m maxic_ws13_hr_hr
## 6 profc_ws11_10m planc_ws7_10m
## k 5
## 1 crosc_ws13_hr_hr
## 2 maxic_ws29_hr_hr
## 3 maxic_ws9_hr_hr
## 4 planc_ws9_hr_hr
## 5 Tangential_Curvature_hr
## 6 slope_ws29_hr_hr
## [1] "OBB error with all predictors of localterrain is 0.509259259259259"
## MeanDecreaseGini MeanDecreaseAccuracy
## Slope_50m 1.1373329 0.006772203
## maxic_ws9_hr_hr 0.9402189 0.005657172
## Slope_10m 0.9382462 0.005540153
## slope_ws11_10m 0.6850574 0.004272931
## slope_DTM_50m_avg_ws9_50m 1.1119233 0.004260792
## slope_ws15_10m 0.7460380 0.004107919
## maxic_ws11_hr_hr 0.7214319 0.003908498
## MaximalCurvature_10m 0.9016103 0.003164278
## slope_ws5_10m 0.7458977 0.003064262
## GeneralCurvature_10m 0.6611461 0.002844019
## [1] "OBB error with all predictors of localterrain is 0.481481481481481"
## [1] "OOB-error: 0.481481481481481 for predictors Slope_50m AND GeneralCurvature_10m"
##
## preds 2 3 4 5
## 2 33 0 0 0
## 3 0 23 0 0
## 4 0 0 35 0
## 5 0 0 0 17
## [1] "Kappa overall = 1"
## [1] "Tau overall = 1"
## [1] "mean quality = 1"
## [1] "The quality is 1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 1"
## [1] "classification error rate with altdata: 0"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of localterrain is 0.472222222222222"
## [1] "OOB-error: 0.425925925925926 for predictors Slope_50m AND GeneralCurvature_10m AND Flow_Line_Curvature_hr AND LongitudinalCurvature_10m AND Minimal_Curvature_hr AND longc_ws5_10m AND Plan_Curvature_hr AND crosc_ws13_hr_hr AND maxic_ws13_hr_hr AND slope_DTM_50m_avg_ws3_50m AND slope_ws29_hr_hr AND Profile_Curvature_50m AND Total_Curvature_50m AND Convexity_50m AND maxic_ws11_10m AND slope_DTM_50m_avg_ws3_50m AND profc_ws15_10m AND profc_ws11_10m AND maxic_ws9_hr_hr AND slope_DTM_50m_avg_ws5_50m AND maxic_ws5_hr_hr AND maxic_ws19_hr_hr AND maxic_ws13_hr_hr AND planc_ws7_10m AND crosc_ws13_hr_hr AND maxic_ws29_hr_hr AND maxic_ws9_hr_hr AND planc_ws9_hr_hr AND Tangential_Curvature_hr AND slope_ws29_hr_hr"
##
## preds 2 3 4 5
## 2 33 0 0 0
## 3 0 23 0 0
## 4 0 0 35 0
## 5 0 0 0 17
## [1] "Kappa overall = 1"
## [1] "Tau overall = 1"
## [1] "mean quality = 1"
## [1] "The quality is 1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 1"
## [1] "classification error rate with altdata: 0"
## [1] "1a.3 - Habitat for soil organisms"
## [1] "Prediction error at end is: 0.675324675324675"
## [2] "Prediction error at end is: 0.694372294372294"
## [3] "Prediction error at end is: 0.648917748917749"
## [4] "Prediction error at end is: 0.602597402597403"
## [5] "Prediction error at end is: 0.528571428571429"
## [6] "Prediction error at end is: 0.567099567099567"
## k 1 k 2
## 1 minic_ws9_hr_hr minic_ws9_hr_hr
## 2 slope_ws19_hr_hr profc_DTM_50m_avg_ws3_50m
## 3 minic_DTM_50m_avg_ws7_50m crosc_ws13_hr_hr
## 4 slope_DTM_50m_avg_ws9_50m dtm_hr_CONVEX_r30
## 5 minic_ws5_hr_hr planc_ws15_hr_hr
## 6 planc_ws9_hr_hr crosc_DTM_50m_avg_ws11_50m
## k 3 k 4
## 1 crosc_DTM_50m_avg_ws7_50m FlowLineCurvature_10m
## 2 crosc_ws15_hr_hr Total_Curvature_50m
## 3 slope_DTM_50m_avg_ws3_50m profc_ws15_10m
## 4 Convexity_50m Convergence_Index_hr
## 5 planc_DTM_50m_avg_ws5_50m slope_DTM_50m_avg_ws7_50m
## 6 CrossSectionalCurvature_50m maxic_ws3_hr_hr
## k 5
## 1 minic_ws9_hr_hr
## 2 slope_ws15_10m
## 3 dtm_hr_CONVEX_r30
## 4 profc_DTM_50m_avg_ws11_50m
## 5 minic_ws23_hr_hr
## 6 crosc_ws15_10m
## [1] "OBB error with all predictors of localterrain is 0.453703703703704"
## MeanDecreaseGini MeanDecreaseAccuracy
## TotalCurvature_10m 1.6823580 0.007306444
## slope_DTM_50m_avg_ws3_50m 1.2725379 0.005536008
## minic_ws5_hr_hr 1.1382578 0.004968556
## maxic_ws15_hr_hr 0.8943811 0.004615553
## maxic_ws23_hr_hr 0.7709114 0.004033272
## Slope_50m 0.8929971 0.003954678
## maxic_DTM_50m_avg_ws9_50m 0.8513740 0.003878230
## slope_ws15_10m 0.9317333 0.003659854
## slope_DTM_50m_avg_ws5_50m 0.9155822 0.003635030
## maxic_DTM_50m_avg_ws5_50m 0.6806564 0.003433799
## [1] "OBB error with all predictors of localterrain is 0.444444444444444"
## [1] "OOB-error: 0.564814814814815 for predictors minic_ws9_hr_hr AND slope_ws19_hr_hr"
##
## preds 1 2 3 4
## 1 24 0 0 0
## 2 0 34 0 0
## 3 0 0 26 0
## 4 0 0 0 24
## [1] "Kappa overall = 1"
## [1] "Tau overall = 1"
## [1] "mean quality = 1"
## [1] "The quality is 1"
## [1] "######### Cramer's V = 1"
## [1] "classification error rate with altdata: 0"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of localterrain is 0.490740740740741"
## [1] "OOB-error: 0.444444444444444 for predictors minic_ws9_hr_hr AND slope_ws19_hr_hr AND minic_DTM_50m_avg_ws7_50m AND slope_DTM_50m_avg_ws9_50m AND minic_ws5_hr_hr AND planc_ws9_hr_hr AND minic_ws9_hr_hr AND profc_DTM_50m_avg_ws3_50m AND crosc_ws13_hr_hr AND dtm_hr_CONVEX_r30 AND planc_ws15_hr_hr AND crosc_DTM_50m_avg_ws11_50m AND crosc_DTM_50m_avg_ws7_50m AND crosc_ws15_hr_hr AND slope_DTM_50m_avg_ws3_50m AND Convexity_50m AND planc_DTM_50m_avg_ws5_50m AND CrossSectionalCurvature_50m AND FlowLineCurvature_10m AND Total_Curvature_50m AND profc_ws15_10m AND Convergence_Index_hr AND slope_DTM_50m_avg_ws7_50m AND maxic_ws3_hr_hr AND minic_ws9_hr_hr AND slope_ws15_10m AND dtm_hr_CONVEX_r30 AND profc_DTM_50m_avg_ws11_50m AND minic_ws23_hr_hr AND crosc_ws15_10m"
##
## preds 1 2 3 4
## 1 24 0 0 0
## 2 0 34 0 0
## 3 0 0 26 0
## 4 0 0 0 24
## [1] "Kappa overall = 1"
## [1] "Tau overall = 1"
## [1] "mean quality = 1"
## [1] "The quality is 1"
## [1] "######### Cramer's V = 1"
## [1] "classification error rate with altdata: 0"
## [1] "1a.4 - Habitat for crops"
## [1] "Prediction error at end is: 0.267099567099567"
## [2] "Prediction error at end is: 0.183982683982684"
## [3] "Prediction error at end is: 0.174458874458874"
## [4] "Prediction error at end is: 0.212121212121212"
## [5] "Prediction error at end is: 0.212121212121212"
## [6] "Prediction error at end is: 0.211688311688312"
## k 1 k 2 k 3
## 1 slope_ws5_hr_hr slope_ws5_10m slope_ws3_hr_hr
## 2 crosc_ws13_hr_hr maxic_ws13_hr_hr longc_ws15_hr_hr
## 3 Slope_hr slope_ws3_10m minic_DTM_50m_avg_ws9_50m
## 4 minic_ws23_hr_hr minic_ws11_10m longc_ws7_10m
## 5 minic_ws5_10m Tangential_Curvature_hr longc_ws3_hr_hr
## 6 planc_ws19_hr_hr maxic_ws7_hr_hr profc_DTM_50m_avg_ws7_50m
## k 4 k 5
## 1 slope_ws3_hr_hr Slope_hr
## 2 minic_ws15_hr_hr longc_ws15_hr_hr
## 3 minic_DTM_50m_avg_ws9_50m slope_ws5_hr_hr
## 4 Convexity_10m profc_ws15_hr_hr
## 5 longc_ws7_hr_hr slope_DTM_50m_avg_ws9_50m
## 6 minic_ws23_hr_hr ProfileCurvature_10m
## [1] "OBB error with all predictors of localterrain is 0.157407407407407"
## MeanDecreaseGini MeanDecreaseAccuracy
## slope_ws3_hr_hr 3.090598 0.02574590
## Slope_hr 3.030718 0.02508142
## slope_ws5_hr_hr 2.709475 0.02110797
## slope_ws23_hr_hr 2.400080 0.01759921
## Slope_10m 2.016602 0.01630050
## slope_ws5_10m 1.773646 0.01471760
## slope_ws19_hr_hr 1.977858 0.01424233
## slope_ws13_hr_hr 2.081435 0.01402077
## slope_ws11_hr_hr 2.180392 0.01332207
## slope_ws7_10m 1.696628 0.01207609
## [1] "OBB error with all predictors of localterrain is 0.185185185185185"
## [1] "OOB-error: 0.12037037037037 for predictors slope_ws5_hr_hr AND crosc_ws13_hr_hr"
##
## preds 3 4 5
## 3 8 0 0
## 4 0 62 0
## 5 0 0 38
## [1] "Kappa overall = 1"
## [1] "Tau overall = 1"
## [1] "mean quality = 1"
## [1] "The quality is 1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 1"
## [1] "classification error rate with altdata: 0"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of localterrain is 0.175925925925926"
## [1] "OOB-error: 0.12037037037037 for predictors slope_ws5_hr_hr AND crosc_ws13_hr_hr AND Slope_hr AND minic_ws23_hr_hr AND minic_ws5_10m AND planc_ws19_hr_hr AND slope_ws5_10m AND maxic_ws13_hr_hr AND slope_ws3_10m AND minic_ws11_10m AND Tangential_Curvature_hr AND maxic_ws7_hr_hr AND slope_ws3_hr_hr AND longc_ws15_hr_hr AND minic_DTM_50m_avg_ws9_50m AND longc_ws7_10m AND longc_ws3_hr_hr AND profc_DTM_50m_avg_ws7_50m AND slope_ws3_hr_hr AND minic_ws15_hr_hr AND minic_DTM_50m_avg_ws9_50m AND Convexity_10m AND longc_ws7_hr_hr AND minic_ws23_hr_hr AND Slope_hr AND longc_ws15_hr_hr AND slope_ws5_hr_hr AND profc_ws15_hr_hr AND slope_DTM_50m_avg_ws9_50m AND ProfileCurvature_10m"
##
## preds 3 4 5
## 3 8 0 0
## 4 0 62 0
## 5 0 0 38
## [1] "Kappa overall = 1"
## [1] "Tau overall = 1"
## [1] "mean quality = 1"
## [1] "The quality is 1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 1"
## [1] "classification error rate with altdata: 0"
## [1] "1c.1 - Average precipitation retention capacity"
## [1] "Prediction error at end is: 0.751082251082251"
## [2] "Prediction error at end is: 0.674458874458874"
## [3] "Prediction error at end is: 0.629437229437229"
## [4] "Prediction error at end is: 0.583982683982684"
## [5] "Prediction error at end is: 0.601298701298701"
## [6] "Prediction error at end is: 0.592207792207792"
## k 1 k 2
## 1 crosc_ws23_hr_hr crosc_DTM_50m_avg_ws9_50m
## 2 planc_ws5_10m profc_ws7_10m
## 3 maxic_DTM_50m_avg_ws3_50m crosc_ws5_10m
## 4 slope_DTM_50m_avg_ws11_50m Convergence_Index_hr
## 5 slope_DTM_50m_avg_ws5_50m planc_DTM_50m_avg_ws9_50m
## 6 dtm_hr_CONVEX_r10 slope_ws15_10m
## k 3 k 4 k 5
## 1 slope_ws29_hr_hr crosc_ws23_hr_hr slope_ws3_hr_hr
## 2 slope_ws5_hr_hr slope_ws11_hr_hr maxic_ws3_10m
## 3 Minimal_Curvature_hr maxic_ws3_10m slope_DTM_50m_avg_ws9_50m
## 4 crosc_DTM_50m_avg_ws11_50m maxic_ws23_hr_hr maxic_ws5_hr_hr
## 5 minic_ws15_10m planc_ws19_hr_hr ProfileCurvature_10m
## 6 TotalCurvature_10m Slope_50m slope_ws15_10m
## [1] "OBB error with all predictors of localterrain is 0.648148148148148"
## MeanDecreaseGini MeanDecreaseAccuracy
## MaximalCurvature_10m 0.7707003 0.002423896
## profc_ws3_hr_hr 0.7067442 0.002401000
## Minimal_Curvature_hr 0.5479792 0.002055238
## Slope_hr 0.6140058 0.002037640
## slope_DTM_50m_avg_ws11_50m 0.7824028 0.001692910
## Flow_Line_Curvature_50m 0.7049653 0.001623029
## longc_ws3_hr_hr 0.7170861 0.001527864
## ProfileCurvature_10m 0.6099676 0.001442390
## maxic_ws5_10m 0.5446579 0.001433610
## minic_ws5_hr_hr 0.8674175 0.001372780
## [1] "OBB error with all predictors of localterrain is 0.62037037037037"
## [1] "OOB-error: 0.546296296296296 for predictors crosc_ws23_hr_hr AND planc_ws5_10m"
##
## preds 1 2 3 4 5
## 1 26 0 0 0 0
## 2 0 45 0 0 0
## 3 0 0 12 0 0
## 4 0 0 0 18 0
## 5 0 0 0 0 7
## [1] "Kappa overall = 1"
## [1] "Tau overall = 1"
## [1] "mean quality = 1"
## [1] "The quality is 1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 1"
## [1] "classification error rate with altdata: 0"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of localterrain is 0.62037037037037"
## [1] "OOB-error: 0.574074074074074 for predictors crosc_ws23_hr_hr AND planc_ws5_10m AND maxic_DTM_50m_avg_ws3_50m AND slope_DTM_50m_avg_ws11_50m AND slope_DTM_50m_avg_ws5_50m AND dtm_hr_CONVEX_r10 AND crosc_DTM_50m_avg_ws9_50m AND profc_ws7_10m AND crosc_ws5_10m AND Convergence_Index_hr AND planc_DTM_50m_avg_ws9_50m AND slope_ws15_10m AND slope_ws29_hr_hr AND slope_ws5_hr_hr AND Minimal_Curvature_hr AND crosc_DTM_50m_avg_ws11_50m AND minic_ws15_10m AND TotalCurvature_10m AND crosc_ws23_hr_hr AND slope_ws11_hr_hr AND maxic_ws3_10m AND maxic_ws23_hr_hr AND planc_ws19_hr_hr AND Slope_50m AND slope_ws3_hr_hr AND maxic_ws3_10m AND slope_DTM_50m_avg_ws9_50m AND maxic_ws5_hr_hr AND ProfileCurvature_10m AND slope_ws15_10m"
##
## preds 1 2 3 4 5
## 1 26 0 0 0 0
## 2 0 45 0 0 0
## 3 0 0 12 0 0
## 4 0 0 0 18 0
## 5 0 0 0 0 7
## [1] "Kappa overall = 1"
## [1] "Tau overall = 1"
## [1] "mean quality = 1"
## [1] "The quality is 1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 1"
## [1] "classification error rate with altdata: 0"
## [1] "1c.1 - Minimum precipitation retention capacity"
## [1] "Prediction error at end is: 0.814718614718615"
## [2] "Prediction error at end is: 0.693073593073593"
## [3] "Prediction error at end is: 0.75021645021645"
## [4] "Prediction error at end is: 0.694372294372294"
## [5] "Prediction error at end is: 0.675324675324675"
## [6] "Prediction error at end is: 0.702597402597403"
## k 1 k 2 k 3
## 1 profc_ws7_10m longc_ws23_hr_hr profc_ws29_hr_hr
## 2 crosc_ws7_hr_hr crosc_ws29_hr_hr maxic_ws15_hr_hr
## 3 planc_ws29_hr_hr minic_ws15_hr_hr minic_ws7_10m
## 4 Total_Curvature_hr crosc_ws5_10m profc_ws3_10m
## 5 crosc_ws3_10m Flow_Line_Curvature_50m minic_ws13_hr_hr
## 6 minic_ws15_hr_hr profc_ws11_10m maxic_ws29_hr_hr
## k 4 k 5
## 1 maxic_ws7_10m planc_ws29_hr_hr
## 2 planc_ws13_hr_hr longc_ws5_10m
## 3 crosc_ws5_10m crosc_ws15_hr_hr
## 4 longc_DTM_50m_avg_ws9_50m crosc_ws15_10m
## 5 CrossSectionalCurvature_hr planc_DTM_50m_avg_ws11_50m
## 6 Tangential_Curvature_50m profc_ws5_10m
## [1] "OBB error with all predictors of localterrain is 0.666666666666667"
## MeanDecreaseGini MeanDecreaseAccuracy
## TotalCurvature_10m 1.6915471 0.004384354
## crosc_ws5_10m 0.6672548 0.003256847
## maxic_ws13_hr_hr 0.6509108 0.002767572
## Total_Curvature_50m 1.2497922 0.002696361
## crosc_ws11_10m 0.6882563 0.002473567
## crosc_ws7_10m 0.4615653 0.002408718
## maxic_ws15_hr_hr 0.7767013 0.001615967
## FlowLineCurvature_10m 0.6296937 0.001611732
## maxic_ws11_hr_hr 0.6001174 0.001608941
## planc_ws11_10m 0.5128324 0.001540164
## [1] "OBB error with all predictors of localterrain is 0.62962962962963"
## [1] "OOB-error: 0.657407407407407 for predictors profc_ws7_10m AND crosc_ws7_hr_hr"
##
## preds 1 2 3 4 5
## 1 12 0 0 0 0
## 2 0 26 0 0 0
## 3 0 0 22 0 0
## 4 0 0 0 26 0
## 5 0 0 0 0 22
## [1] "Kappa overall = 1"
## [1] "Tau overall = 1"
## [1] "mean quality = 1"
## [1] "The quality is 1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 1"
## [1] "classification error rate with altdata: 0"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of localterrain is 0.666666666666667"
## [1] "OOB-error: 0.490740740740741 for predictors profc_ws7_10m AND crosc_ws7_hr_hr AND planc_ws29_hr_hr AND Total_Curvature_hr AND crosc_ws3_10m AND minic_ws15_hr_hr AND longc_ws23_hr_hr AND crosc_ws29_hr_hr AND minic_ws15_hr_hr AND crosc_ws5_10m AND Flow_Line_Curvature_50m AND profc_ws11_10m AND profc_ws29_hr_hr AND maxic_ws15_hr_hr AND minic_ws7_10m AND profc_ws3_10m AND minic_ws13_hr_hr AND maxic_ws29_hr_hr AND maxic_ws7_10m AND planc_ws13_hr_hr AND crosc_ws5_10m AND longc_DTM_50m_avg_ws9_50m AND CrossSectionalCurvature_hr AND Tangential_Curvature_50m AND planc_ws29_hr_hr AND longc_ws5_10m AND crosc_ws15_hr_hr AND crosc_ws15_10m AND planc_DTM_50m_avg_ws11_50m AND profc_ws5_10m"
##
## preds 1 2 3 4 5
## 1 12 0 0 0 0
## 2 0 26 0 0 0
## 3 0 0 22 0 0
## 4 0 0 0 26 0
## 5 0 0 0 0 22
## [1] "Kappa overall = 1"
## [1] "Tau overall = 1"
## [1] "mean quality = 1"
## [1] "The quality is 1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 1"
## [1] "classification error rate with altdata: 0"
## [1] "1c.2 - Retention capacity for heavy precipitation events"
## [1] "Prediction error at end is: 0.471428571428571"
## [2] "Prediction error at end is: 0.352380952380952"
## [3] "Prediction error at end is: 0.287878787878788"
## [4] "Prediction error at end is: 0.296536796536797"
## [5] "Prediction error at end is: 0.296536796536797"
## [6] "Prediction error at end is: 0.296969696969697"
## k 1 k 2
## 1 Minimal_Curvature_50m maxic_ws23_hr_hr
## 2 Longitudinal_Curvature_hr minic_ws29_hr_hr
## 3 longc_DTM_50m_avg_ws3_50m Convexity_50m
## 4 MaximalCurvature_10m slope_DTM_50m_avg_ws9_50m
## 5 maxic_ws13_hr_hr Plan_Curvature_hr
## 6 minic_ws15_10m MaximalCurvature_10m
## k 3 k 4
## 1 crosc_ws7_10m Convexity_50m
## 2 slope_DTM_50m_avg_ws9_50m longc_ws15_10m
## 3 crosc_ws5_hr_hr FlowLineCurvature_10m
## 4 minic_ws5_hr_hr minic_ws3_hr_hr
## 5 profc_DTM_50m_avg_ws5_50m profc_ws11_10m
## 6 Plan_Curvature_50m Flow_Line_Curvature_hr
## k 5
## 1 maxic_ws23_hr_hr
## 2 minic_ws23_hr_hr
## 3 Minimal_Curvature_50m
## 4 maxic_DTM_50m_avg_ws11_50m
## 5 profc_DTM_50m_avg_ws11_50m
## 6 minic_ws15_10m
## [1] "OBB error with all predictors of localterrain is 0.277777777777778"
## MeanDecreaseGini MeanDecreaseAccuracy
## FlowLineCurvature_10m 0.9386443 0.006205879
## crosc_DTM_50m_avg_ws5_50m 0.8284482 0.005603521
## minic_ws7_10m 0.7772594 0.005294088
## minic_ws29_hr_hr 0.8524506 0.004889761
## crosc_ws11_10m 0.6672203 0.003382759
## maxic_ws7_10m 0.3935556 0.003287878
## maxic_ws23_hr_hr 0.3907311 0.003249979
## minic_ws23_hr_hr 0.6889049 0.003240223
## minic_ws19_hr_hr 0.5945955 0.003059828
## longc_ws5_10m 0.5116574 0.002863669
## [1] "OBB error with all predictors of localterrain is 0.305555555555556"
## [1] "OOB-error: 0.240740740740741 for predictors Minimal_Curvature_50m AND Longitudinal_Curvature_hr"
##
## preds 1 2 3 4 5
## 1 75 0 0 0 0
## 2 0 11 0 0 0
## 3 0 0 6 0 0
## 4 0 0 0 8 0
## 5 0 0 0 0 8
## [1] "Kappa overall = 1"
## [1] "Tau overall = 1"
## [1] "mean quality = 1"
## [1] "The quality is 1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 1"
## [1] "classification error rate with altdata: 0"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of localterrain is 0.277777777777778"
## [1] "OOB-error: 0.25 for predictors Minimal_Curvature_50m AND Longitudinal_Curvature_hr AND longc_DTM_50m_avg_ws3_50m AND MaximalCurvature_10m AND maxic_ws13_hr_hr AND minic_ws15_10m AND maxic_ws23_hr_hr AND minic_ws29_hr_hr AND Convexity_50m AND slope_DTM_50m_avg_ws9_50m AND Plan_Curvature_hr AND MaximalCurvature_10m AND crosc_ws7_10m AND slope_DTM_50m_avg_ws9_50m AND crosc_ws5_hr_hr AND minic_ws5_hr_hr AND profc_DTM_50m_avg_ws5_50m AND Plan_Curvature_50m AND Convexity_50m AND longc_ws15_10m AND FlowLineCurvature_10m AND minic_ws3_hr_hr AND profc_ws11_10m AND Flow_Line_Curvature_hr AND maxic_ws23_hr_hr AND minic_ws23_hr_hr AND Minimal_Curvature_50m AND maxic_DTM_50m_avg_ws11_50m AND profc_DTM_50m_avg_ws11_50m AND minic_ws15_10m"
##
## preds 1 2 3 4 5
## 1 75 0 0 0 0
## 2 0 11 0 0 0
## 3 0 0 6 0 0
## 4 0 0 0 8 0
## 5 0 0 0 0 8
## [1] "Kappa overall = 1"
## [1] "Tau overall = 1"
## [1] "mean quality = 1"
## [1] "The quality is 1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 1"
## [1] "classification error rate with altdata: 0"
## [1] "1c.3 - groundwater reformation rate"
## [1] "Prediction error at end is: 0.788311688311688"
## [2] "Prediction error at end is: 0.73030303030303"
## [3] "Prediction error at end is: 0.693506493506494"
## [4] "Prediction error at end is: 0.700865800865801"
## [5] "Prediction error at end is: 0.674025974025974"
## [6] "Prediction error at end is: 0.701731601731602"
## k 1 k 2
## 1 slope_ws15_hr_hr LongitudinalCurvature_10m
## 2 profc_ws23_hr_hr maxic_ws7_10m
## 3 maxic_ws19_hr_hr maxic_ws3_hr_hr
## 4 slope_DTM_50m_avg_ws5_50m maxic_ws29_hr_hr
## 5 longc_DTM_50m_avg_ws3_50m slope_DTM_50m_avg_ws5_50m
## 6 minic_ws29_hr_hr GeneralCurvature_10m
## k 3 k 4
## 1 minic_ws15_10m minic_ws11_hr_hr
## 2 FlowLineCurvature_10m maxic_DTM_50m_avg_ws5_50m
## 3 dtm_hr_CONVEX_r5 crosc_ws13_hr_hr
## 4 planc_ws11_10m slope_DTM_50m_avg_ws9_50m
## 5 minic_DTM_50m_avg_ws3_50m planc_ws3_10m
## 6 crosc_ws7_10m maxic_ws9_hr_hr
## k 5
## 1 minic_ws15_10m
## 2 maxic_ws19_hr_hr
## 3 longc_ws3_hr_hr
## 4 maxic_ws15_hr_hr
## 5 maxic_DTM_50m_avg_ws5_50m
## 6 planc_ws7_10m
## [1] "OBB error with all predictors of localterrain is 0.611111111111111"
## MeanDecreaseGini MeanDecreaseAccuracy
## slope_ws5_hr_hr 0.9091755 0.004505613
## slope_ws23_hr_hr 0.8713569 0.003693752
## slope_ws19_hr_hr 0.6994681 0.003065583
## Profile_Curvature_hr 0.8664606 0.003013326
## slope_ws7_10m 0.7242359 0.002625295
## slope_ws15_hr_hr 0.6548519 0.002503940
## minic_ws3_hr_hr 0.9761192 0.002479712
## slope_ws5_10m 0.8099342 0.002418275
## Slope_10m 0.7466590 0.002312560
## Minimal_Curvature_hr 0.5855733 0.002036061
## [1] "OBB error with all predictors of localterrain is 0.611111111111111"
## [1] "OOB-error: 0.62037037037037 for predictors slope_ws15_hr_hr AND profc_ws23_hr_hr"
##
## preds 1 2 3 4 5
## 1 14 0 0 0 0
## 2 0 39 0 0 0
## 3 0 0 14 0 0
## 4 0 0 0 19 0
## 5 0 0 0 0 22
## [1] "Kappa overall = 1"
## [1] "Tau overall = 1"
## [1] "mean quality = 1"
## [1] "The quality is 1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 1"
## [1] "classification error rate with altdata: 0"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of localterrain is 0.611111111111111"
## [1] "OOB-error: 0.564814814814815 for predictors slope_ws15_hr_hr AND profc_ws23_hr_hr AND maxic_ws19_hr_hr AND slope_DTM_50m_avg_ws5_50m AND longc_DTM_50m_avg_ws3_50m AND minic_ws29_hr_hr AND LongitudinalCurvature_10m AND maxic_ws7_10m AND maxic_ws3_hr_hr AND maxic_ws29_hr_hr AND slope_DTM_50m_avg_ws5_50m AND GeneralCurvature_10m AND minic_ws15_10m AND FlowLineCurvature_10m AND dtm_hr_CONVEX_r5 AND planc_ws11_10m AND minic_DTM_50m_avg_ws3_50m AND crosc_ws7_10m AND minic_ws11_hr_hr AND maxic_DTM_50m_avg_ws5_50m AND crosc_ws13_hr_hr AND slope_DTM_50m_avg_ws9_50m AND planc_ws3_10m AND maxic_ws9_hr_hr AND minic_ws15_10m AND maxic_ws19_hr_hr AND longc_ws3_hr_hr AND maxic_ws15_hr_hr AND maxic_DTM_50m_avg_ws5_50m AND planc_ws7_10m"
##
## preds 1 2 3 4 5
## 1 14 0 0 0 0
## 2 0 39 0 0 0
## 3 0 0 14 0 0
## 4 0 0 0 19 0
## 5 0 0 0 0 22
## [1] "Kappa overall = 1"
## [1] "Tau overall = 1"
## [1] "mean quality = 1"
## [1] "The quality is 1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 1"
## [1] "classification error rate with altdata: 0"
## [1] "1c.4 - Potential for providing nutrients for plants"
## [1] "Prediction error at end is: 0.529004329004329"
## [2] "Prediction error at end is: 0.371428571428571"
## [3] "Prediction error at end is: 0.325108225108225"
## [4] "Prediction error at end is: 0.324675324675325"
## [5] "Prediction error at end is: 0.306493506493506"
## [6] "Prediction error at end is: 0.288311688311688"
## k 1 k 2
## 1 minic_ws5_hr_hr crosc_ws5_hr_hr
## 2 longc_DTM_50m_avg_ws5_50m profc_DTM_50m_avg_ws11_50m
## 3 maxic_ws15_hr_hr planc_DTM_50m_avg_ws11_50m
## 4 profc_ws29_hr_hr minic_ws19_hr_hr
## 5 maxic_ws15_10m Plan_Curvature_50m
## 6 maxic_ws11_10m minic_ws7_hr_hr
## k 3 k 4
## 1 longc_ws29_hr_hr TangentialCurvature_10m
## 2 minic_ws7_10m slope_ws11_10m
## 3 longc_ws7_10m minic_DTM_50m_avg_ws7_50m
## 4 crosc_ws5_hr_hr slope_DTM_50m_avg_ws7_50m
## 5 profc_ws29_hr_hr CrossSectionalCurvature_50m
## 6 CrossSectionalCurvature_50m maxic_DTM_50m_avg_ws5_50m
## k 5
## 1 longc_ws19_hr_hr
## 2 Total_Curvature_50m
## 3 minic_DTM_50m_avg_ws7_50m
## 4 dtm_hr_CONVEX_r30
## 5 planc_DTM_50m_avg_ws11_50m
## 6 minic_ws15_hr_hr
## [1] "OBB error with all predictors of localterrain is 0.305555555555556"
## MeanDecreaseGini MeanDecreaseAccuracy
## minic_ws11_hr_hr 1.3896970 0.011297624
## minic_ws9_hr_hr 1.3926995 0.007597971
## crosc_ws5_hr_hr 1.1285306 0.007037121
## minic_ws7_hr_hr 1.5603885 0.006705559
## TotalCurvature_10m 1.2511444 0.006160002
## minic_ws5_hr_hr 0.9093794 0.005050313
## planc_ws5_hr_hr 0.6506936 0.003306029
## TangentialCurvature_10m 0.3848811 0.003082114
## MaximalCurvature_10m 0.4665715 0.003028937
## maxic_ws7_10m 0.7210332 0.003016545
## [1] "OBB error with all predictors of localterrain is 0.305555555555556"
## [1] "OOB-error: 0.296296296296296 for predictors minic_ws5_hr_hr AND longc_DTM_50m_avg_ws5_50m"
##
## preds 1 3 5
## 1 65 0 0
## 3 0 11 0
## 5 0 0 32
## [1] "Kappa overall = 1"
## [1] "Tau overall = 1"
## [1] "mean quality = 1"
## [1] "The quality is 1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 1"
## [1] "classification error rate with altdata: 0"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of localterrain is 0.305555555555556"
## [1] "OOB-error: 0.25 for predictors minic_ws5_hr_hr AND longc_DTM_50m_avg_ws5_50m AND maxic_ws15_hr_hr AND profc_ws29_hr_hr AND maxic_ws15_10m AND maxic_ws11_10m AND crosc_ws5_hr_hr AND profc_DTM_50m_avg_ws11_50m AND planc_DTM_50m_avg_ws11_50m AND minic_ws19_hr_hr AND Plan_Curvature_50m AND minic_ws7_hr_hr AND longc_ws29_hr_hr AND minic_ws7_10m AND longc_ws7_10m AND crosc_ws5_hr_hr AND profc_ws29_hr_hr AND CrossSectionalCurvature_50m AND TangentialCurvature_10m AND slope_ws11_10m AND minic_DTM_50m_avg_ws7_50m AND slope_DTM_50m_avg_ws7_50m AND CrossSectionalCurvature_50m AND maxic_DTM_50m_avg_ws5_50m AND longc_ws19_hr_hr AND Total_Curvature_50m AND minic_DTM_50m_avg_ws7_50m AND dtm_hr_CONVEX_r30 AND planc_DTM_50m_avg_ws11_50m AND minic_ws15_hr_hr"
##
## preds 1 3 5
## 1 65 0 0
## 3 0 11 0
## 5 0 0 32
## [1] "Kappa overall = 1"
## [1] "Tau overall = 1"
## [1] "mean quality = 1"
## [1] "The quality is 1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 1"
## [1] "classification error rate with altdata: 0"
## [1] "1c.5 - Potential as a CO2 sink"
## [1] "Prediction error at end is: 0.522943722943723"
## [2] "Prediction error at end is: 0.486580086580087"
## [3] "Prediction error at end is: 0.421212121212121"
## [4] "Prediction error at end is: 0.42987012987013"
## [5] "Prediction error at end is: 0.401731601731602"
## [6] "Prediction error at end is: 0.354112554112554"
## k 1 k 2 k 3
## 1 slope_ws19_hr_hr slope_ws19_hr_hr CrossSectionalCurvature_hr
## 2 PlanCurvature_10m minic_DTM_50m_avg_ws11_50m TotalCurvature_10m
## 3 crosc_ws5_hr_hr planc_ws5_10m FlowLineCurvature_10m
## 4 maxic_ws3_hr_hr planc_ws13_hr_hr planc_DTM_50m_avg_ws7_50m
## 5 TotalCurvature_10m Tangential_Curvature_50m Total_Curvature_50m
## 6 planc_ws3_10m crosc_ws13_hr_hr profc_DTM_50m_avg_ws5_50m
## k 4 k 5
## 1 crosc_ws5_hr_hr crosc_ws19_hr_hr
## 2 slope_ws13_hr_hr profc_ws7_10m
## 3 minic_ws15_hr_hr Tangential_Curvature_50m
## 4 longc_ws7_10m maxic_ws15_10m
## 5 slope_ws15_hr_hr maxic_ws29_hr_hr
## 6 PlanCurvature_10m TangentialCurvature_10m
## [1] "OBB error with all predictors of localterrain is 0.351851851851852"
## MeanDecreaseGini MeanDecreaseAccuracy
## TotalCurvature_10m 1.8085695 0.011117328
## maxic_ws15_hr_hr 1.6261493 0.009747365
## minic_ws11_hr_hr 0.9883106 0.006768397
## MaximalCurvature_10m 1.0611950 0.006696107
## minic_ws3_hr_hr 0.9256518 0.005123103
## minic_ws5_hr_hr 1.1148332 0.004956734
## CrossSectionalCurvature_hr 0.7767361 0.004492110
## slope_ws5_10m 0.6749511 0.004464976
## maxic_ws23_hr_hr 0.6333793 0.004412928
## Total_Curvature_hr 1.0246174 0.004146658
## [1] "OBB error with all predictors of localterrain is 0.37037037037037"
## [1] "OOB-error: 0.37037037037037 for predictors slope_ws19_hr_hr AND PlanCurvature_10m"
##
## preds 1 2 3 4 5
## 1 51 0 0 0 0
## 2 0 1 0 0 0
## 3 0 0 18 0 0
## 4 0 0 0 33 0
## 5 0 0 0 0 5
## [1] "Kappa overall = 1"
## [1] "Tau overall = 1"
## [1] "mean quality = 1"
## [1] "The quality is 1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 1"
## [1] "classification error rate with altdata: 0"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of localterrain is 0.324074074074074"
## [1] "OOB-error: 0.314814814814815 for predictors slope_ws19_hr_hr AND PlanCurvature_10m AND crosc_ws5_hr_hr AND maxic_ws3_hr_hr AND TotalCurvature_10m AND planc_ws3_10m AND slope_ws19_hr_hr AND minic_DTM_50m_avg_ws11_50m AND planc_ws5_10m AND planc_ws13_hr_hr AND Tangential_Curvature_50m AND crosc_ws13_hr_hr AND CrossSectionalCurvature_hr AND TotalCurvature_10m AND FlowLineCurvature_10m AND planc_DTM_50m_avg_ws7_50m AND Total_Curvature_50m AND profc_DTM_50m_avg_ws5_50m AND crosc_ws5_hr_hr AND slope_ws13_hr_hr AND minic_ws15_hr_hr AND longc_ws7_10m AND slope_ws15_hr_hr AND PlanCurvature_10m AND crosc_ws19_hr_hr AND profc_ws7_10m AND Tangential_Curvature_50m AND maxic_ws15_10m AND maxic_ws29_hr_hr AND TangentialCurvature_10m"
##
## preds 1 2 3 4 5
## 1 51 0 0 0 0
## 2 0 1 0 0 0
## 3 0 0 18 0 0
## 4 0 0 0 33 0
## 5 0 0 0 0 5
## [1] "Kappa overall = 1"
## [1] "Tau overall = 1"
## [1] "mean quality = 1"
## [1] "The quality is 1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 1"
## [1] "classification error rate with altdata: 0"
## [1] "1d.1 - Potential for retention of heavy metals"
## [1] "Prediction error at end is: 0.638095238095238"
## [2] "Prediction error at end is: 0.5004329004329"
## [3] "Prediction error at end is: 0.435930735930736"
## [4] "Prediction error at end is: 0.427705627705628"
## [5] "Prediction error at end is: 0.408658008658009"
## [6] "Prediction error at end is: 0.399134199134199"
## k 1 k 2
## 1 slope_ws15_10m longc_ws23_hr_hr
## 2 minic_ws7_10m slope_ws29_hr_hr
## 3 dtm_hr_CONVEX_r30 maxic_ws7_10m
## 4 longc_DTM_50m_avg_ws5_50m CrossSectionalCurvature_hr
## 5 planc_DTM_50m_avg_ws5_50m crosc_ws5_hr_hr
## 6 crosc_ws15_hr_hr crosc_ws13_hr_hr
## k 3 k 4
## 1 minic_ws29_hr_hr Total_Curvature_50m
## 2 Flow_Line_Curvature_50m CrossSectionalCurvature_hr
## 3 maxic_ws5_hr_hr minic_ws3_hr_hr
## 4 planc_ws3_hr_hr dtm_hr_CONVEX_r30
## 5 dtm_hr_CONVEX_r20 crosc_DTM_50m_avg_ws7_50m
## 6 FlowLineCurvature_10m TotalCurvature_10m
## k 5
## 1 slope_DTM_50m_avg_ws3_50m
## 2 maxic_ws29_hr_hr
## 3 slope_ws23_hr_hr
## 4 Total_Curvature_hr
## 5 crosc_ws13_hr_hr
## 6 Minimal_Curvature_hr
## [1] "OBB error with all predictors of localterrain is 0.361111111111111"
## MeanDecreaseGini MeanDecreaseAccuracy
## TotalCurvature_10m 1.4197148 0.007170865
## slope_ws3_10m 0.9638191 0.005567992
## minic_ws5_hr_hr 1.3393282 0.005490887
## maxic_ws15_hr_hr 0.9104475 0.005090625
## Slope_10m 1.1650392 0.005079061
## slope_ws5_10m 0.8172995 0.004812107
## Slope_50m 0.6501166 0.004447252
## slope_ws13_hr_hr 0.8421275 0.004091610
## slope_ws23_hr_hr 0.8940395 0.004091067
## crosc_ws13_hr_hr 0.6049590 0.003628726
## [1] "OBB error with all predictors of localterrain is 0.37037037037037"
## [1] "OOB-error: 0.37962962962963 for predictors slope_ws15_10m AND minic_ws7_10m"
##
## preds 1 2 3 4 5
## 1 42 0 0 0 0
## 2 0 7 0 0 0
## 3 0 0 10 0 0
## 4 0 0 0 7 0
## 5 0 0 0 0 42
## [1] "Kappa overall = 1"
## [1] "Tau overall = 1"
## [1] "mean quality = 1"
## [1] "The quality is 1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 1"
## [1] "classification error rate with altdata: 0"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of localterrain is 0.37037037037037"
## [1] "OOB-error: 0.333333333333333 for predictors slope_ws15_10m AND minic_ws7_10m AND dtm_hr_CONVEX_r30 AND longc_DTM_50m_avg_ws5_50m AND planc_DTM_50m_avg_ws5_50m AND crosc_ws15_hr_hr AND longc_ws23_hr_hr AND slope_ws29_hr_hr AND maxic_ws7_10m AND CrossSectionalCurvature_hr AND crosc_ws5_hr_hr AND crosc_ws13_hr_hr AND minic_ws29_hr_hr AND Flow_Line_Curvature_50m AND maxic_ws5_hr_hr AND planc_ws3_hr_hr AND dtm_hr_CONVEX_r20 AND FlowLineCurvature_10m AND Total_Curvature_50m AND CrossSectionalCurvature_hr AND minic_ws3_hr_hr AND dtm_hr_CONVEX_r30 AND crosc_DTM_50m_avg_ws7_50m AND TotalCurvature_10m AND slope_DTM_50m_avg_ws3_50m AND maxic_ws29_hr_hr AND slope_ws23_hr_hr AND Total_Curvature_hr AND crosc_ws13_hr_hr AND Minimal_Curvature_hr"
##
## preds 1 2 3 4 5
## 1 42 0 0 0 0
## 2 0 7 0 0 0
## 3 0 0 10 0 0
## 4 0 0 0 7 0
## 5 0 0 0 0 42
## [1] "Kappa overall = 1"
## [1] "Tau overall = 1"
## [1] "mean quality = 1"
## [1] "The quality is 1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 1"
## [1] "classification error rate with altdata: 0"
## [1] "1d.2 - Potential for transforming organic contaminants"
## [1] "Prediction error at end is: 0.657142857142857"
## [2] "Prediction error at end is: 0.647186147186147"
## [3] "Prediction error at end is: 0.593073593073593"
## [4] "Prediction error at end is: 0.62034632034632"
## [5] "Prediction error at end is: 0.603463203463203"
## [6] "Prediction error at end is: 0.601298701298701"
## k 1 k 2 k 3
## 1 crosc_ws19_hr_hr MaximalCurvature_10m crosc_ws15_hr_hr
## 2 Minimal_Curvature_hr maxic_ws3_hr_hr crosc_ws19_hr_hr
## 3 maxic_ws29_hr_hr MinimalCurvature_10m maxic_ws15_hr_hr
## 4 minic_DTM_50m_avg_ws9_50m longc_ws15_10m minic_DTM_50m_avg_ws5_50m
## 5 crosc_ws3_hr_hr dtm_hr_CONVEX_r5 longc_DTM_50m_avg_ws9_50m
## 6 slope_DTM_50m_avg_ws3_50m planc_ws15_hr_hr Plan_Curvature_hr
## k 4 k 5
## 1 maxic_ws7_10m crosc_ws13_hr_hr
## 2 maxic_ws23_hr_hr crosc_DTM_50m_avg_ws3_50m
## 3 longc_ws11_10m planc_DTM_50m_avg_ws5_50m
## 4 profc_ws7_hr_hr crosc_ws3_10m
## 5 minic_ws13_hr_hr CrossSectionalCurvature_50m
## 6 longc_ws7_10m minic_ws5_hr_hr
## [1] "OBB error with all predictors of localterrain is 0.638888888888889"
## MeanDecreaseGini MeanDecreaseAccuracy
## TotalCurvature_10m 1.3895773 0.005646626
## MaximalCurvature_10m 1.3814653 0.005120462
## minic_ws5_hr_hr 0.8674868 0.003467020
## slope_ws19_hr_hr 0.7628082 0.003461461
## slope_ws29_hr_hr 0.6657954 0.003205129
## slope_ws5_10m 0.7759054 0.002965819
## minic_ws7_hr_hr 0.8871311 0.002909758
## slope_ws7_10m 0.8566747 0.002885358
## maxic_ws9_hr_hr 0.7843232 0.002630408
## maxic_ws15_hr_hr 0.9504995 0.002571349
## [1] "OBB error with all predictors of localterrain is 0.601851851851852"
## [1] "OOB-error: 0.555555555555556 for predictors crosc_ws19_hr_hr AND Minimal_Curvature_hr"
##
## preds 1 2 3 4 5
## 1 10 0 0 0 0
## 2 0 22 0 0 0
## 3 0 0 28 0 0
## 4 0 0 0 9 0
## 5 0 0 0 0 39
## [1] "Kappa overall = 1"
## [1] "Tau overall = 1"
## [1] "mean quality = 1"
## [1] "The quality is 1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 1"
## [1] "classification error rate with altdata: 0"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of localterrain is 0.601851851851852"
## [1] "OOB-error: 0.564814814814815 for predictors crosc_ws19_hr_hr AND Minimal_Curvature_hr AND maxic_ws29_hr_hr AND minic_DTM_50m_avg_ws9_50m AND crosc_ws3_hr_hr AND slope_DTM_50m_avg_ws3_50m AND MaximalCurvature_10m AND maxic_ws3_hr_hr AND MinimalCurvature_10m AND longc_ws15_10m AND dtm_hr_CONVEX_r5 AND planc_ws15_hr_hr AND crosc_ws15_hr_hr AND crosc_ws19_hr_hr AND maxic_ws15_hr_hr AND minic_DTM_50m_avg_ws5_50m AND longc_DTM_50m_avg_ws9_50m AND Plan_Curvature_hr AND maxic_ws7_10m AND maxic_ws23_hr_hr AND longc_ws11_10m AND profc_ws7_hr_hr AND minic_ws13_hr_hr AND longc_ws7_10m AND crosc_ws13_hr_hr AND crosc_DTM_50m_avg_ws3_50m AND planc_DTM_50m_avg_ws5_50m AND crosc_ws3_10m AND CrossSectionalCurvature_50m AND minic_ws5_hr_hr"
##
## preds 1 2 3 4 5
## 1 10 0 0 0 0
## 2 0 22 0 0 0
## 3 0 0 28 0 0
## 4 0 0 0 9 0
## 5 0 0 0 0 39
## [1] "Kappa overall = 1"
## [1] "Tau overall = 1"
## [1] "mean quality = 1"
## [1] "The quality is 1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 1"
## [1] "classification error rate with altdata: 0"
## [1] "1d.3 - Potential as filter and buffer for organic contaminants"
## [1] "Prediction error at end is: 0.0554112554112554"
## [2] "Prediction error at end is: 0.0463203463203463"
## [3] "Prediction error at end is: 0.0372294372294372"
## [4] "Prediction error at end is: 0.0372294372294372"
## [5] "Prediction error at end is: 0.0372294372294372"
## [6] "Prediction error at end is: 0.0372294372294372"
## k 1 k 2
## 1 Convexity_10m PlanCurvature_10m
## 2 planc_ws5_10m planc_ws9_hr_hr
## 3 Convexity_50m slope_ws13_hr_hr
## 4 Maximal_Curvature_50m maxic_DTM_50m_avg_ws9_50m
## 5 Minimal_Curvature_50m slope_DTM_50m_avg_ws3_50m
## 6 Plan_Curvature_50m slope_ws15_10m
## k 3 k 4
## 1 CrossSectionalCurvature_50m longc_ws11_10m
## 2 PlanCurvature_10m longc_ws7_10m
## 3 slope_ws13_hr_hr profc_ws5_10m
## 4 Flow_Line_Curvature_50m Minimal_Curvature_50m
## 5 maxic_ws7_hr_hr crosc_ws7_hr_hr
## 6 planc_ws11_hr_hr maxic_ws23_hr_hr
## k 5
## 1 Convexity_10m
## 2 planc_ws9_hr_hr
## 3 PlanCurvature_10m
## 4 minic_ws29_hr_hr
## 5 CrossSectionalCurvature_50m
## 6 minic_ws23_hr_hr
## [1] "OBB error with all predictors of localterrain is 0.037037037037037"
## MeanDecreaseGini MeanDecreaseAccuracy
## profc_DTM_50m_avg_ws5_50m 0.05348294 0.0015212168
## planc_DTM_50m_avg_ws5_50m 0.14353614 0.0014809356
## longc_ws15_10m 0.05926461 0.0013900200
## planc_DTM_50m_avg_ws11_50m 0.21200915 0.0012811719
## profc_DTM_50m_avg_ws9_50m 0.06087598 0.0011735909
## profc_ws11_10m 0.05079702 0.0011057148
## profc_ws15_10m 0.03807431 0.0010371631
## minic_ws29_hr_hr 0.07496966 0.0009558680
## profc_DTM_50m_avg_ws7_50m 0.10393737 0.0008745406
## longc_DTM_50m_avg_ws7_50m 0.05407145 0.0008603118
## [1] "OBB error with all predictors of localterrain is 0.037037037037037"
## [1] "OOB-error: 0.037037037037037 for predictors Convexity_10m AND planc_ws5_10m"
##
## preds 4 5
## 4 4 0
## 5 0 104
## [1] "Kappa overall = 1"
## [1] "Tau overall = 1"
## [1] "mean quality = 1"
## [1] "The quality is 1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 0.870192307692308"
## [1] "classification error rate with altdata: 0"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of localterrain is 0.037037037037037"
## [1] "OOB-error: 0.0185185185185185 for predictors Convexity_10m AND planc_ws5_10m AND Convexity_50m AND Maximal_Curvature_50m AND Minimal_Curvature_50m AND Plan_Curvature_50m AND PlanCurvature_10m AND planc_ws9_hr_hr AND slope_ws13_hr_hr AND maxic_DTM_50m_avg_ws9_50m AND slope_DTM_50m_avg_ws3_50m AND slope_ws15_10m AND CrossSectionalCurvature_50m AND PlanCurvature_10m AND slope_ws13_hr_hr AND Flow_Line_Curvature_50m AND maxic_ws7_hr_hr AND planc_ws11_hr_hr AND longc_ws11_10m AND longc_ws7_10m AND profc_ws5_10m AND Minimal_Curvature_50m AND crosc_ws7_hr_hr AND maxic_ws23_hr_hr AND Convexity_10m AND planc_ws9_hr_hr AND PlanCurvature_10m AND minic_ws29_hr_hr AND CrossSectionalCurvature_50m AND minic_ws23_hr_hr"
##
## preds 4 5
## 4 4 0
## 5 0 104
## [1] "Kappa overall = 1"
## [1] "Tau overall = 1"
## [1] "mean quality = 1"
## [1] "The quality is 1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 0.870192307692308"
## [1] "classification error rate with altdata: 0"
## [1] "1d.4 - Potential for retention of water-soluble contaminants"
## [1] "Prediction error at end is: 0.795238095238095"
## [2] "Prediction error at end is: 0.600865800865801"
## [3] "Prediction error at end is: 0.572727272727273"
## [4] "Prediction error at end is: 0.554545454545455"
## [5] "Prediction error at end is: 0.526406926406926"
## [6] "Prediction error at end is: 0.535930735930736"
## k 1 k 2
## 1 Minimal_Curvature_hr crosc_ws19_hr_hr
## 2 crosc_ws3_10m CrossSectionalCurvature_hr
## 3 crosc_ws7_hr_hr crosc_ws23_hr_hr
## 4 Slope_50m profc_DTM_50m_avg_ws9_50m
## 5 planc_DTM_50m_avg_ws3_50m crosc_DTM_50m_avg_ws3_50m
## 6 Profile_Curvature_50m profc_ws5_hr_hr
## k 3 k 4
## 1 GeneralCurvature_10m maxic_ws13_hr_hr
## 2 FlowLineCurvature_10m slope_ws29_hr_hr
## 3 crosc_ws3_hr_hr Tangential_Curvature_hr
## 4 crosc_ws15_hr_hr LongitudinalCurvature_10m
## 5 longc_ws3_10m crosc_ws11_hr_hr
## 6 MinimalCurvature_10m longc_DTM_50m_avg_ws3_50m
## k 5
## 1 crosc_ws13_hr_hr
## 2 MinimalCurvature_10m
## 3 crosc_DTM_50m_avg_ws9_50m
## 4 crosc_ws15_hr_hr
## 5 crosc_DTM_50m_avg_ws11_50m
## 6 profc_ws7_hr_hr
## [1] "OBB error with all predictors of localterrain is 0.5"
## MeanDecreaseGini MeanDecreaseAccuracy
## slope_ws11_hr_hr 0.8151606 0.004513121
## slope_ws3_10m 0.8861537 0.003939477
## TotalCurvature_10m 0.9456823 0.003589144
## slope_ws19_hr_hr 0.6407110 0.003355717
## slope_ws23_hr_hr 0.6931134 0.003346353
## slope_DTM_50m_avg_ws5_50m 0.7497832 0.003238108
## Slope_10m 0.7486609 0.003200858
## crosc_ws7_10m 0.7558984 0.003080825
## crosc_ws3_10m 0.5450793 0.002890011
## Minimal_Curvature_hr 0.8398992 0.002520341
## [1] "OBB error with all predictors of localterrain is 0.5"
## [1] "OOB-error: 0.537037037037037 for predictors Minimal_Curvature_hr AND crosc_ws3_10m"
##
## preds 1 2 3 4 5
## 1 47 0 0 0 0
## 2 0 14 0 0 0
## 3 0 0 17 0 0
## 4 0 0 0 22 0
## 5 0 0 0 0 8
## [1] "Kappa overall = 1"
## [1] "Tau overall = 1"
## [1] "mean quality = 1"
## [1] "The quality is 1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 1"
## [1] "classification error rate with altdata: 0"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of localterrain is 0.490740740740741"
## [1] "OOB-error: 0.481481481481481 for predictors Minimal_Curvature_hr AND crosc_ws3_10m AND crosc_ws7_hr_hr AND Slope_50m AND planc_DTM_50m_avg_ws3_50m AND Profile_Curvature_50m AND crosc_ws19_hr_hr AND CrossSectionalCurvature_hr AND crosc_ws23_hr_hr AND profc_DTM_50m_avg_ws9_50m AND crosc_DTM_50m_avg_ws3_50m AND profc_ws5_hr_hr AND GeneralCurvature_10m AND FlowLineCurvature_10m AND crosc_ws3_hr_hr AND crosc_ws15_hr_hr AND longc_ws3_10m AND MinimalCurvature_10m AND maxic_ws13_hr_hr AND slope_ws29_hr_hr AND Tangential_Curvature_hr AND LongitudinalCurvature_10m AND crosc_ws11_hr_hr AND longc_DTM_50m_avg_ws3_50m AND crosc_ws13_hr_hr AND MinimalCurvature_10m AND crosc_DTM_50m_avg_ws9_50m AND crosc_ws15_hr_hr AND crosc_DTM_50m_avg_ws11_50m AND profc_ws7_hr_hr"
##
## preds 1 2 3 4 5
## 1 47 0 0 0 0
## 2 0 14 0 0 0
## 3 0 0 17 0 0
## 4 0 0 0 22 0
## 5 0 0 0 0 8
## [1] "Kappa overall = 1"
## [1] "Tau overall = 1"
## [1] "mean quality = 1"
## [1] "The quality is 1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 1"
## [1] "classification error rate with altdata: 0"
## [1] "1d.5 - Potential as buffer for acidic contaminants"
## [1] "Prediction error at end is: 0.759307359307359"
## [2] "Prediction error at end is: 0.684848484848485"
## [3] "Prediction error at end is: 0.674458874458874"
## [4] "Prediction error at end is: 0.674891774891775"
## [5] "Prediction error at end is: 0.637662337662338"
## [6] "Prediction error at end is: 0.647186147186147"
## k 1 k 2 k 3
## 1 crosc_ws29_hr_hr crosc_ws23_hr_hr planc_ws7_hr_hr
## 2 FlowLineCurvature_10m planc_ws11_10m maxic_ws23_hr_hr
## 3 planc_DTM_50m_avg_ws11_50m longc_ws15_10m planc_DTM_50m_avg_ws3_50m
## 4 crosc_ws3_hr_hr planc_ws3_10m crosc_ws19_hr_hr
## 5 crosc_ws13_hr_hr profc_ws15_10m crosc_ws5_hr_hr
## 6 planc_ws11_10m crosc_ws15_10m planc_ws23_hr_hr
## k 4 k 5
## 1 minic_DTM_50m_avg_ws9_50m slope_ws3_10m
## 2 crosc_ws5_hr_hr planc_ws5_hr_hr
## 3 slope_DTM_50m_avg_ws9_50m crosc_ws19_hr_hr
## 4 Slope_50m crosc_ws5_hr_hr
## 5 longc_ws5_hr_hr maxic_ws11_10m
## 6 longc_ws15_hr_hr slope_ws13_hr_hr
## [1] "OBB error with all predictors of localterrain is 0.611111111111111"
## MeanDecreaseGini MeanDecreaseAccuracy
## maxic_DTM_50m_avg_ws11_50m 1.3473532 0.004675391
## longc_DTM_50m_avg_ws11_50m 0.8828896 0.003151859
## Total_Curvature_50m 0.9422908 0.003136906
## minic_ws11_hr_hr 0.8122509 0.002817321
## crosc_ws5_hr_hr 0.9757649 0.002705295
## Slope_10m 0.8251192 0.002654011
## slope_ws3_10m 0.6112843 0.002107568
## planc_ws23_hr_hr 0.6790395 0.002102952
## slope_ws15_hr_hr 0.6160903 0.002085516
## minic_ws5_10m 0.5756878 0.002005237
## [1] "OBB error with all predictors of localterrain is 0.611111111111111"
## [1] "OOB-error: 0.638888888888889 for predictors crosc_ws29_hr_hr AND FlowLineCurvature_10m"
##
## preds 1 2 3 4 5
## 1 39 0 0 0 0
## 2 0 24 0 0 0
## 3 0 0 20 0 0
## 4 0 0 0 16 0
## 5 0 0 0 0 9
## [1] "Kappa overall = 1"
## [1] "Tau overall = 1"
## [1] "mean quality = 1"
## [1] "The quality is 1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 1"
## [1] "classification error rate with altdata: 0"
## [1] "###################### WITH PREDICTORS from the FW SELECTION ###################"
## [1] "OBB error with all predictors of localterrain is 0.592592592592593"
## [1] "OOB-error: 0.555555555555556 for predictors crosc_ws29_hr_hr AND FlowLineCurvature_10m AND planc_DTM_50m_avg_ws11_50m AND crosc_ws3_hr_hr AND crosc_ws13_hr_hr AND planc_ws11_10m AND crosc_ws23_hr_hr AND planc_ws11_10m AND longc_ws15_10m AND planc_ws3_10m AND profc_ws15_10m AND crosc_ws15_10m AND planc_ws7_hr_hr AND maxic_ws23_hr_hr AND planc_DTM_50m_avg_ws3_50m AND crosc_ws19_hr_hr AND crosc_ws5_hr_hr AND planc_ws23_hr_hr AND minic_DTM_50m_avg_ws9_50m AND crosc_ws5_hr_hr AND slope_DTM_50m_avg_ws9_50m AND Slope_50m AND longc_ws5_hr_hr AND longc_ws15_hr_hr AND slope_ws3_10m AND planc_ws5_hr_hr AND crosc_ws19_hr_hr AND crosc_ws5_hr_hr AND maxic_ws11_10m AND slope_ws13_hr_hr"
##
## preds 1 2 3 4 5
## 1 39 0 0 0 0
## 2 0 24 0 0 0
## 3 0 0 20 0 0
## 4 0 0 0 16 0
## 5 0 0 0 0 9
## [1] "Kappa overall = 1"
## [1] "Tau overall = 1"
## [1] "mean quality = 1"
## [1] "The quality is 1"
## Warning in chisq.test(CM): Chi-squared approximation may be incorrect
## [1] "######### Cramer's V = 1"
## [1] "classification error rate with altdata: 0"